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Abstract. We study Darboux-type transformations associated with the fo- 
cusing nonlinear Schrodinger equation (NLS_) and their effect on spectral 
properties of the underlying Lax operator. The latter is a formally ^-self- 
adjoint (but non-self-adjoint) Dirac-typc differential expression of the form 



satisfying JM{q)J = M(q)*, where J is denned by J = (?J)C, and C 
denotes the antilinear conjugation map in C 2 , C(a,b) T = (a,b) T , a,b G C. As 
one of our principal results we prove that under the most general hypothesis 
q £ (R) on q, the maximally defined operator D(q) generated by M(q) 
is actually j7-self-adjoint in L 2 (K) 2 . Moreover, we establish the existence 
of Weyl-Titchmarsh-type solutions •) € L 2 ([R, oo)) 2 and *_(*;,•) G 

L 2 ((— oo, R]) for all R e R of M(q)9±(z) = z*±(z) for z in the resolvent set 
of D. 

The Darboux transformations considered in this paper are the analog of 
the double commutation procedure familiar in the KdV and Schrodinger op- 
erator contexts. As in the corresponding case of Schrodinger operators, the 
Darboux transformations in question guarantee that the resulting potentials 
q are locally nonsingular. Moreover, we prove that the construction of N- 
soliton NLS- potentials q^> with respect to a general NLS_ background 
potential q £ Lj- (R), associated with the Dirac-type operators D(q( N ^) and 
D(q), respectively, amounts to the insertion of TV complex conjugate pairs of 
L 2 (M) 2 -eigenvalues {z±,zi,...,zn,~zn} into the spectrum a{D(q)) of D(q), 
leaving the rest of the spectrum (especially, the essential spectrum tr c (D(g))) 
invariant, that is, 



These results are obtained by establishing the existence of bounded transfor- 
mation operators which intertwine the background Dirac operator D(q) and 
the Dirac operator D(q( N ^) obtained after TV Darboux-type transformations. 



Date: October, 2002. 

1991 Mathematics Subject Classification. Primary: 34L05, 34L40, 35Q51, 35Q55. Secondary: 
34B20, 47A10. 

Key words and phrases. Dirac operator, focusing nonlinear Schrodinger equation, transforma- 
tion operators, Darboux transformations, J-self-adjointness, Weyl-Titchmarsh solutions. 




dx 



a{D{q { - N ^)) = a(D(q)) U {z u Si, . . .,z N ,z N }, 
<T (D(qW)) =o- c (D(q)). 
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1. Introduction 



Various methods of inserting eigenvalues in spectral gaps of one-dimensional 
Schrodinger operators H(q) associated with differential expressions of the type 



in L 2 (M) (or in L 2 ((a, oo)), a > — oo), with q real-valued and locally integrable, have 
attracted an enormous amount of attention. This is due to their prominent role 
in diverse fields such as the inverse scattering approach, supersymmetric quantum 
mechanics, level comparison theorems, as a tool to construct soliton solutions of the 
Korteweg-de Vries (KdV) hierarchy relative to (general) KdV background solutions, 
and in connection with Backlund transformations for the KdV hierarchy. The 
literature on this subject is too extensive to go into details here, but we refer to 
the detailed accounts given in [14], [15], [16, App. G], [17], [18] and the references 
cited therein. Historically, these methods of inserting eigenvalues go back to Jacobi 
[24] and Darboux [9] with decisive later contributions by Crum [8] , Schmincke [40] , 
and, especially, Deift [10]. 

Two particular methods turned out to be of special importance: The single 
commutation method, also called the Crum-Darboux method [8] , [9] (actually going 
back at least to Jacobi [24]) and the double commutation method, to be found, for 
instance, in the seminal work of Gel'fand and Levitan [13]. (The latter can be 
obtained by a composition of two separate single commutation steps, explaining 
the name double commutation.) 

The single commutation method, although very simply implemented, has the 
distinct disadvantage of relying on positivity properties of certain solutions ip of 
H(q)ip — \il>, which confines its applicability to the insertion of eigenvalues below 
the spectrum of H{q) (assuming H{q) to be bounded from below). A complete 
spectral characterization of this method has been provided by Deift [10] (see also 
[40]) on the basis of unitary equivalence of A* A\ ker ^± and AA* |k er (A*)- L f° r a 
densely defined closed linear operator A in a (complex, separable) Hilbcrt space. 

The double commutation method on the other hand, allows one to insert eigen- 
values into any spectral gap of H(q). Although relatively simply implemented 
also, a complete spectral characterization of the double commutation method for 
Schrodinger-type operators was more recently achieved in [14] on the basis of Weyl- 
Titchmarsh m-function techniques and subsequently in [18] (for general Sturm- 
Liouville operators on arbitrary intervals) using a functional analytic approach 
based on the notion of (intertwining) transformation operators. 

In this paper we concentrate on the analog of the double commutation method for 
Dirac-type operators associated with the Lax operator for the focusing nonlinear 
Schrodinger (NLS_) hierarchy Assuming q to be locally integrable, the Dirac- 
type operator corresponding to the Lax differential expression in the NLS_ case is 
associated with the 2x2 matrix-valued differential expression 



for x E K (cf., e.g., [12, Part I, Sect. 1.2], [46], and [47]). The maximally defined 
Dirac-type operator associated with M(q) in the (two-component) Hilbert space 
L 2 (R) 2 will then be denoted by D(q). By way of contrast, the corresponding 




(1.2) 
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(1.3) 



(formally self-adjoint) Lax differential expression for the defocusing NLS+ case is 
given by 

'£ ~i 

As it turns out there is no natural analog of the single commutation method 
for the Dirac operators associated with the focusing and defocusing nonlinear 
Schrodinger hierarchies (NLS±). However, the complexified version of the NLS± 
hierarchies, the Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy, supports two nat- 
ural analogs of the single commutation method. In order to briefly describe them, 
we recall that the Dirac-type Lax differential expression associated with the AKNS 
hierarchy is given by 

M(p,q) = i(*> ~l) (1.4) 

V V dx/ 

(cf. e.g., [1] and [16, Ch. 3]) in terms of two locally integrable coefficients p,q on 
R. The focusing (NLS_) and defocusing (NLS + ) nonlinear Schrodinger hierarchies 
are then associated with the constraints 



NLS±: p(x) = ±q(x), (1.5) 

respectively. In this paper we will concentrate on the focusing NLS_ case only. 

The two analogs of the single commutation method for the AKNS case, which 
are usually called elementary Darboux transformations, can then be described as 
follows. Suppose 

M(p,q)^(z u x) = z^izux), ^{z u x) = {ip 1 {z 1 ,x),ip 2 (z 1 ,x)) T , (z u x) eCxl, 

(1.6) 

and 

M{jp,qj®{z u x) = Z!^(zi,x), V(z u x) = {■tp 1 {z 1 ,x),i> 2 (zi,x)) T , {z u x) eCxl. 

(1.7) 

Then the two elementary Darboux transformations in the AKNS context are given 
by (cf. [27], [28]) 

(p,q) i-» {p Zl ,q Zl ), (1.8) 

where 



p Zl {x) = -1p2{zi,x)/lpi(zi,x), 

q Zl {x) = q'(x) - ip 2 (zi,x)/i> 1 (z 1 ,x)q(x)' 2 + 2iz 1 q(x), 



(1.9) 



and 
where 



(p,q) ^ {p Zl ,q Zl ), (1.10) 

(1.11) 



p Zl {x) = -p'{x) + 4>i(zi, x) /ip2(zi, x)p{xf + 2izip(x), 

q Zl {x) = ^l{zx,x)/^2{zi,x). 

Similar to the case of Schrodinger operators, the analog of the double commuta- 
tion method for Dirac-type operators associated with (1.4) is then obtained by an 
appropriate composition of the two elementary Darboux transformations (1.9) and 
(1.11). This two-step procedure is denoted by 

(P,9)"(P£*,<&) (1-12) 
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and leads to (cf., e.g., [27], [28], [34, Sect. 4.2], [39]) 

P { zlzM =P( X ) -2i(«i -zi)M*i> x)ij 2 (zi,x)/W(^( Zl , x), a;)), i 

(Here W%F, G) denotes the Wronskian of F, G £ C 2 ). In contrast to (1.8), (1.9) 
and (1.10), (1.11), the two-step procedure (1.12), (1.13) with z\ = z\, is compat- 
ible with the NLS± cases and one explicitly obtains the following Darboux-type 
transformation in the NLS- case, 



The transformation (1.14) for Dirac-type operators assoiated with (1.2) in the NLS_ 
context represents the analog of the double commutation method for Schrodingcr 
operators and leads to locally nonsingular NLS_ potentials qi\\ assuming q to be 
free of local singularities. 

By analogy to the KdV and Schrodinger operator case, one expects the NLS_ 
potential q z \ (x) to produce an eigenvalue at the spectral point z\ for the associ- 
ated Dirac operator .D(gi^), assuming z\ to be a point in the resolvent set of the 
"background" operator D(q). Actually, by a simple symmetry consideration, one 
expects a pair of eigenvalues (z 1; zT) in the point spectrum of D(qi]^). To prove 
this fact and to show that the remaining spectral characteristics (especially, the 
essential spectrum of D(q)) remain invariant under the Darboux-type transforma- 
tion (1.14), is the principal purpose of this paper. More precisely, if we denote by 
qi^}..,z N the NLS- potential obtained after an A-fold iteration of the Darboux-type 
transformation and by D(q z ^'.. tZN ) the resulting Dirac operator, we will prove that 

°( D (<l£}..,z„)) = °( D (<l)) U ■ ■ ■ . Z N ,Z N }, (1.15) 

<r p (D(qi»l iZN )) = a p (D(q)) U { Zl , Zl , . . . , z Nl z N }, (1.16) 

<y C {D(q^l ZN ))=a c (D(q)). (1.17) 

(Here <j(S), <j p (S), and a e (S) denote the spectrum, point spectrum, and essential 
spectrum of a densely defined closed operator 5 in a complex separable Hilbcrt 
space H, cf. Section 5 for more details on spectra, etc.) Actually, we will go a 
step beyond (1.15)— (1.17) and establish the existence of bounded transformation 
operators which intertwine D(q z ^}..^ ZN ) and D{q). 

When trying to embark on proving results of the type (1.15)— (1. 17) for Dirac-type 
operators associated with the differential expression (1.2) in the NLS_ context, one 
finds oneself at a distinct disadvantage when compared to the case of Schrodingcr 
operators with real-valued potentials: While L in (1.1) is formally self-adjoint for 
q real-valued, M(q) in (1.2) is never self-adjoint (except, in the trivial case q = 0). 
As a consequence, the original approach to a complete spectral characterization of 
the double commutation method for Schrodinger operators in [14], based on Weyl- 
Titchmarsh theory and hence on spectral theory, is doomed from the start as there 
simply is no spectral theory for non-self-adjoint Dirac-type operators associated 
with (1.2) under our general hypothesis q s L\ oc {M^). That leaves one with only 
one possible line of attack, the analog of the transformation operator approach 
developed in the general Sturm-Liouville context in [18]. As it will turn out in 
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Section 6, this approach is indeed successful although it requires more sophisticated 
and elaborate arguments compared to those in [18]. 

While the differential expression M{q) in (1.2) is never formally self-adjoint (if 
q ^ 0), it is, however, formally ^-self-adjoint, that is, 

JM(q)J = M(q)*. (1.18) 

Here J is defined by 

J=Q J)c, (1.19) 

and C denotes the antilinear conjugation map in C 2 , 

C(a,b) T = (a,&) T , a,6eC. (1.20) 

As one of our principal results in this paper we will prove in Section 3 that under 
the most general hypothesis q E L 1 1 oc (R), the maximally defined Dirac operator 
D(q) associated with M(q) is in fact j7-self-adjoint, 

JD(q)J = D(q)* (= D(-q)). (1.21) 

As an aside we should mention that the corresponding maximally defined Lax 
operator associated with the defocusing NLS+ differential expression (1.3) is in fact 
sclf-adjoint assuming q E Lj oc (R) only (this is proved in the references mentioned in 
Section 3). This should be contrasted with the case of one-dimensional Schrodingcr 
differential expressions L(q) in (1.1) (the Lax differential expression associated with 
the KdV hierarchy). If q E L 1 1 oc (K) in (1.1) is real-valued, then L(q) is formally 
sclf-adjoint but the maximally defined operator H(q) in L 2 (M.) associated with L(q) 
may not be self-adjoint. The latter situation occurs precisely when L(q) is in the 
limit circle case (as opposed to the limit point case) at +oo and/or — oo (cf. [7, Ch. 
9]). This is in sharp contrast to the focusing (respectively, defocusing) NLS case 
where D(q) is always ^-self-adjoint (respectively, self-adjoint). 

Summarizing, we derive the following principal new results in this paper, assum- 
ing the optimal condition q E Ll oc (M.): 

— J'-self-adjointness of D(q). 

— The existence of Weyl-Titchmarsh-type solutions of M(q) , i> = for all 
z E p(D(q)). 

— The existence and boundcdness of transformation operators intertwining the 
operators D(q( N )) and D(q). 

— A spectral analysis of NLS_ Darboux transformations (cf. (1.15)-(1.17)). 

Finally we briefly describe the content of each section. Section 2 first introduces 
our main notation and then proceeds to a review of Darboux transformations for 
AKNS and NLS_ systems. Section 3 is devoted to a proof of the j7-sclf-adjointness 
property (1.21) of D(q), assuming q E L 1 1 oc (M) only. Section 4 constructs eigenval- 
ues of D(qi\ ) ) at pairs zi,zi, z\ E p(D(qj) and associated L 2 (K) 2 -eigenfunctions. 
Section 5 derives some basic spectral properties of general Dirac-type operators 
D(q) and establishes the existence of Weyl-Titchmarsh-type solutions associated 
with M(q). This shows a remarkable similarity to sclf-adjoint systems and appears 
to be without precedent in this non-self-adjoint context. Our final Section 6 es- 
tablishes the existence of bounded transformation operators intertwining D(qi\^ 
and D(q) and then employs these transformation operators to prove the spectral 
properties (1.15)-(1.17). 
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All results in the principal part of this paper, Sections 3-6, are proved under the 
optimal condition q <G L 1 1 oc (]R). Moreover, practically all results in Sections 3-6 are 
new as long as one goes beyond bounded or periodic potentials q. In particular, 
Theorem 6.14 (characterizing transformation operators) and Theorem 6.15 (proving 
(1.15)— (1.17)) appear to be the first of their kind under any assumptions on q. 

In this paper we confine ourselves to a stationary (i.e., time-independent) ap- 
proach only. Applications to the time-dependent focusing NLS_ equation and to 
nonlinear optics will be made in a subsequent paper [3] . 



2. Darboux-type Transformations for AKNS and NLS_ Systems 

In this section we take a close look at Darboux-type transformations for non- 
sclf-adjoint Dirac-typc differential expressions M(q) (cf. (2.3)) applicable to AKNS 
systems, with special emphasis on the case of the focusing nonlinear Schrodingcr 
equation NLS_ (cf. (2.7)). 

Throughout this paper we use the following notation: ' = d/dx; for a matrix A 
with complex-valued entries, A T denotes the transposed matrix, A the matrix with 
complex conjugate entries, and A* = A = A T the adjoint matrix. Occasionally, 
we use the following 2x2 matrices, 

/2 =G ?)• cti =C o)' ff3 =G , 2n 

/ A (l 0\ /0 0\ [ ' ! 

a4= {-l oj' a5= {o 0)> a6 = \0 I)' 

If A = (a\, a 2 ) T and B = {bi 1 b 2 ) T arc 2x 1 column- vectors, then (A, B) C 2 = B* A = 
a\b\ + a 2 b 2 denotes the usual scalar product in C 2 , ||^4||c 2 = (l a i| 2 + l a 2| 2 ) 1/ ^ 2 the 
associated norm, and A 1 - = (a 2 , — ffli) = A T the 1x2 row- vector perpendicular 
to A (in the sense that A^A = 0). We also use the notation 

W(A, B) = B^A = -A^B = ai b 2 - a 2 b x = (A,~B^) C2 (2.2) 

for the Wronskian of A and B. The space of 2 x 2 matrices with entries in C 
is denoted by C 2x2 and the operator norm of a 2 x 2 matrix A induced by the 
usual norm in C 2 is denoted by ||A||c2 X 2. In the following, fl C R denotes an 
open subset of M and ACi oc (^) 2xm , m — 1,2, denotes the set of 2 x to matrices 
with locally absolutely continuous entries on Q (for 2x1 columns we use the 
corresponding notation ACi oc (f7) 2 ). We define 1 LP(n) 2xm and Lf oc (f7) 2xm , to = 
1,2, to consist of 2 x to matrices with entries in L P (Q) and L^ oc (Q), respectively. 
In the special case fi = R and F, G G L 2 (M.) 2 , the scalar product of F and G 
is denoted by (F,G)^2 = J R dx (F(x), G(x))c2 with associated norm of F given 

1 /2 

by ||.F||i,2 = (J R dx \\F(x) ||^, 2 ) • Finally, the open complex upper (respectively, 
lower) half-plane is denoted by C+ (respectively, C_); the domain, range, and 
kernel (null space) of a linear operator T are denoted by dom(T), ran(T), and 
ker(T), respectively. 

Hypothesis 2.1. Let S1C1 open and assume p,q e Ll oc (Q). 



For brevity we write L p (fi) for L p (f2;dz) and suppress the Lebesgue measure dx whenever 
possible. 
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Assuming Hypothesis 2.1 and z G C, we introduce the 2x2 matrix U{z,p,q) 
and the 2x2 matrix- valued differential expression M(p, q) by 

U(z,p,q)=(-* Z M(p,q)=i(& _~|) . (2.3) 

The functions p and g in (2.3) are referred to as AKNS potentials due to the fact 
that M(p, q) is the Lax differential expression associated with the AKNS hierarchy 
(see, e.g., [1] and [16, Ch. 3]). The particularly important special case p = —q 
will be referred to as the NLS_ case (due to the obvious connection of (2.3) with 
the zero curvature representation and the Lax operator for the focusing nonlinear 
Schrodinger equation, sec, e.g., [12, Part 1, Sect. 1.2], [46], and [47]), and then q is 
called an NLS_ potential. For given z G C, £1 C R, and AKNS potentials (p, q), a 
function &(z, •) <G AC\ oc (Q) 2 is called a z-wave function associated with (p, q) on 
if ty'(z,x) — U(z,p, q)^(z 7 x) holds for a.e. ieO, that is, if W — (ipi,ip 2 ) T satisfies 
the following first-order system of differential equations 

tp[(z,x) = -izipi(z,x) + q(x)ip 2 (z,x), tp' 2 (z,x) = izip 2 {z,x) + p(x)ipi (z, x) 

(2.4) 

for a.e. x £ ft. Equivalently, ^f(z) — (ipi(z),ip 2 (z)) T G ACi oc (Jl) 2 is a z-wave 
function associated with (p,q) on f2 if and only if M(p, q)^(z) = z^f(z) on £1 in 
the distributional sense. If for some z G C, *(z) and $(z) are z-wave functions 
associated with (p, q) on 0, their Wronskian is well-known to be constant with 
respect to x G fl, 

-^W(^(z,x),^(z,x)) = 0, a: en. (2.5) 

More generally, if ^f(zi) and <E>(z 2 ) are z\- and Z2-wave functions associated with 
(p, q) on Ct, then 

-^W(^(z 1 ,x),^(z 2 ,x)) = i(z 2 - z 1 )[ipi(z 1 ,x)(f> 2 {z 2l x) +ip 2 (z ll x)(j) 1 (z 2 ,x)] 

= i(z 2 - zi)*(zi, x) T cri$(z 2 , x), (2.6) 
f(zi) = (t/>i(zi),V>2(zi)) T , $(z 2 ) = (0i(z 2 ),0 2 (^ 2 )) T , z ly z 2 G C, xeQ. 
In the NLS_ case p = — q we use the notation 

■<**-(4& t% «M* 4) <"> 

instead of (2.3), and we then call any distributional solution ^(z) of M{q)^{z) = 
z^E'(z) an NLS_ z-wave function associated with q. 

If ^(z) is a z-wave function associated with (p,q), and ^(z,x ) = for some 
x G fi, then ^(z,x) vanishes identically for all x in an open neighborhood of x n 
in f2 by the unique solvability of the Cauchy problem for (2.4). Therefore, we will 
always assume that ^(z, x) ^ for all x G fi. 

If q = (and analogously if p = 0) a.e. on fi, then the system (2.4) decomposes 
and yields 

tpi(z, x) = C\ exp(— izx), 

(f x \ ^) 

/ dx' p(x') exp(— 2izx') I exp(izx) + C 2 exp(izx) 

for some constants C\, C 2 G C, xq G O. 
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For further use we collect some simple consequences of (2.4). First we introduce 
the antilinear (i.e., conjugate linear) operator J defined by 

J = a 1 C = (J J)c, J 2 =h, (2.9) 

with C the antilinear conjugation map 

C(a,b) T = (a,b) T , a,beC. (2.10) 
Moreover, we introduce the antilinear involution /C defined by 

JC = <j 4 C=^° 1 J)c, JC 2 = -h. (2.11) 

We also note that 

\\JF\\ C , = \\F\\ C ,, \\KF\\ (? = \\F\\ (? , FeC 2 . (2.12) 

Lemma 2.2. Assume Hypothesis 2.1, z G C, and suppose 'J(z) = (tpi(z), ip2(z)) 
is a z-wave function associated with (p, q) on fl. Then the following assertions hold. 

(i) (f>(z,x) — —ip2(z,x)/ipi(z,x) satisfies the Riccati-type equation 

-(f)' (z, x) + q(x)(f>(z, x) 2 + 2iz(f>(z, x) - p(x) = (2.13) 
on the set tpi(z, x) ^ 0}. 

(ii) <p(z,x) — tpi(z,x)/ijj2(z,x) satisfies the Riccati-type equation 

—<p'(z, x) — p(x)ip(z, x) 2 — 2iz(p(z, x) + q(x) = (2-14) 
on the set {x G O | ip2(z, x) ^ 0}. 

Assume in addition the NLS_ case p = —q. Then the following assertions hold. 
(Hi) If^(z) = (ipi(z),tf>2(z)) T is a z-wave function associated with q then 

K9{z) = a 4 C^(z) = (Ml), -Mz)) T = (*(«) x )* (2.15) 
is a z-wave function associated with q and 

M(q) = -KM{q)K. (2.16) 

(iv) The following identity holds 

(||*(*,aOI&)' = 2Im(z)[|Vi(^^)| 2 - M*,x)\ 2 ]- (2-17) 

(v) M(q) and M(q)* are formally unitarily equivalent in the sense that 

M(q)* = a 3 M(q)a 3 = M(-q). (2.18) 

In addition, M(q) is formally J -self- adjoint in the following sense 

JM(q)J = M(q)*. (2.19) 

Consider AKNS potentials p,q E Ll oc (Q). Fix z,z\,z\ G C, a zi-wave function 
^(zi), and a ii-wave function &(zi) associated with (p,q). Our objective is to 
construct new potentials p^\ q^ E (f^ 1 )) for some open set C f2, and the 
corresponding z-wave functions associated with (p^\q^) on fl^K In the NLS_ 
case p = — q we choose z\ — ~z\ and ^f(zi) — K,^(z\), see Lemma 2.2 (iii). 



SPECTRAL ANALYSIS FOR NLS_ DARBOUX TRANSFORMATIONS 



9 



Remark 2.3. Let T G ACi oc {Q) 2x2 , A, B e L 1 (fl) 2x2 for some open subset Q C R, 
suppose that the following identity holds 

T'{x) + T(x)A(x) - B(x)T(x) = for a.e. i£S], (2.20) 

and assume that $ € ACi oc (^) 2 satisfies the first-order system $' = A<& on fl. 
Then the function <I> (1) , defined by <I> (1) = T$, satisfies G ACi oc (^) 2 and the 
first-order system = £?$W on 0. 

Lemma 2.4. Assume Hypothesis 2.1 and z,Z\ G C. In addition, suppose ^(z\) = 
(ipi(zi) , ip 2 ( z i)) T is a zi-wave function associated with (p,q) on ft and introduce 

n zi ={xen\Mzi,x)^0}. (2.21) 

Define p Zl , q Zl , T (q,p Zl ), and T(z, q,p Zl ) on Cl Zl by 

PzA x ) = -^2{zi,x) /^i{zi,x), q Zl (x) = q'(x) + p Zl (x)q(x) 2 + 2iz 1 q(x), (2.22) 

r ( X ,q,P Zl ) = -U q{X M X) 9 ^), (2-23) 



2 V fciW 1 

t(z,x,q,p Zl ) = i(z - z^os +t (x,q,p Zl ). (2.24) 

Then T = T satisfies (2.20) on Q, Zl with A and B given by 

A(z,x) = U(z,p,q), B(z,x,z-i) = U(z,p Zl ,q Zl ), xeCl Zl . (2.25) 

Proof. Using Lemma 2.2 (i) one verifies that T = T satisfies (2.20) with A = 
U(zi,p, q) and B = U(zi,p Zl ,q Zl ). Since U{z) — U(zi) — i(z — 21)0-3, the conclusion 
follows. □ 

Similarly, by Lemma 2.2 (ii), one has the following analogous result. 

Lemma 2.5. Assume Hypothesis 2.1 and let z,z\ G C. In addition, suppose that 
ty(zi) — (^1(^1), -0 2 (ii)) T is a zi-wave function associated with (p,q) on £1 and 
introduce 

Q Zl ={xen\Mzi,x)^0}. (2.26) 
Definep Zl , q Zl , t (p,q Zl ), and t(z,p, q Zl ) on tt Zl by 
p Zl (x) = -p'(x) +q Zl {x)p(x) 2 +2izip(x), q Zl (x) = tpx (zi, x)/i> 2 (zi, x), (2.27) 

r (x,p,q Zl )= 1 -(-* V (2.28) 



2 \p(x) -p(x)q zi (x) 

T(z,x,p,q Zl ) =i(z- z 1 )a 6 +t (x,p,q Zl ). (2.29) 

Then T = T satisfies (2.20) on Cl Zl with A and B given by 

A{z,x) = U(z,p,q), B(z,x,zi) = U(z,p Zl ,q Zl ), xGQ Zl . (2.30) 

The Darboux-type transformations characterized by (2.24) and (2.29) are also 
called elementary Darboux transformations. They have been discussed, for in- 
stance, in [27] and [28]. In the special context of algebro-geometric AKNS solu- 
tions, the effect of elementary Darboux transformations on the underlying compact 
hyperelliptic curve (in connection with the insertion and deletion of eigenvalues as 
well as the isospectral case) was studied in detail in [15], [16, App. G] (see also [19], 
[20]). 

Next, we will construct the transformation matrix T(z, \l/(zi), \l/(zi)) that sat- 
isfies equation (2.20) with A(z,x) = U(z,p,q) and B(z,x) = U(z,p^ q z ^ Zl ) as 
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the product of f (z,p Zl , q Zl ) and T(z,q,p Zl ). Since we will choose z\ — zj in the 
NLS_ context, we omit the zi-dependence inpW, q^\ c&W, etc., in the NLS_ 
case in the following. 

Theorem 2.6. Assume Hypothesis 2.1. 

(i) Suppose z,z\,z\ G C, and assume that ^(zi) — (tpi(zi), tp2(zi)) T and ^(zi) = 

ip2(zi)) T are z\- and z\-wave functions, respectively, associated with (p,q) 
on ft. In addition, introduce 

^ i ] ii ={xen\W^(z 1 ,x)^(z 1 ,x))^0}. (2.31) 
Define pW Zl , q™ gl , andrf,,*^),*^)) on Sl^ gi by 

p£*» = p(x) - 2i{zt - z 1 )Mzi,x)Mz 1 ,x)/W(^{z 1 ,x)^{S 1 ,x)), (2.32) 
& ( x ) = 9(x) - 2i(zt - z 1 )ij 1 (zi,x)Mz 1 ,x)/W(^{z 1 ,x)^{2 1 ,x)), (2.33) 
r(«,a:,*^i),*(«i)) = - l -zh - ^W^izux),^^^))- 1 

x (^^(ii,^)*^!,^)- 1 - zi^{z 1 ,x)^/{z 1 ,x) 1 -). (2.34) 
Then T satisfies the first-order system 

r'(z, x, *(ii)) + T(z, a;, ^{z{))U{z,p, q) 

~ U{zJ z l]- Zl ,q { Zl ]~ Zl )r(z, x,*(*i), = (2.35) 

a.e. on fi^^ ■ TTjms, ifT(z) is a z -wave function associated with (p,q) on Ct, then 
T z \] Zl ( z )> defined by 

T Z \] Z1 (*> = T(z, *(*i))T(s, a;), (2.36) 

is a z -wave function associated with (p^ Zl , Q z ^ Zl ) on ^ Z i Zl ■ 

(ii) Assume the NLS_ case p = —q and z,z\ € C. Then, for z\ — ~zj, ^f(zi) — 
K?f!(z{) = (ip 2 (zi),-ipi(zi)) , and 

ftW = {xen\ W(*(z!,x), KM(zux)) # 0}, (2.37) 
formulas (2.33)-(2.34) simplify to 

qiWx) = q(x)+4lm(z 1 )^ 1 (z 1 ,x)ip 2 (z ll x)\\^(z 1 ,x)\\-i, (2.38) 
P { ^(x) ^ p(x) - 4lm(z 1 )Mzi,x)Mzi,x)\\^(z 1 ,x)\\^ 



= -q z \\x), (2.39) 
r(z, x, *(*!), /Ctt(*i)) = -(i/2)(2 - Zl )I 2 

+ Im(^)||*(^, x)\\-ilC*(z u x)*{z u x) 1 - (2.40) 

for x £ . In particular, if T(z) is a z-wave function associated with q on Q, 
then Tfi^(z), defined by 

T«(z, x) = T{z, x, ^(z 1 ),JC^>(z 1 ))T(z, x), (2.41) 

is a z-wave function associated with q z ^ on (which may vanish identically 

w.r.t. x e (l^zi , cf. Remark 2.8). 
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Proof. First, we use Lemma 2.4 for z = z\ and ^{z\) to construct p Zl , q Zl , 
f(zi,q,p Zl ) as in (2.22)-(2.24). Define = f (zi, q,p Zl )^{zx). By Lemma 2.4 

and Remark 2.3 we conclude that & Zl (zi) is a zVwave function associated with 
(p Zl ,q Zl ). Moreover, 



^ Zl (Si,x) = i(Si - zi)ipi(zi,x) 

~ \ (Pz 1 (x)lpl(z 1 ,x) +1p2(zi,x)) Q 



(2.42) 



where p Zl is defined in (2.22). We now apply Lemma 2.5 replacing ^(zi) by ^> Zl (zi) 
and (p,q) by (p Zl ,q Zl ). Then q ZlM = ^i, Zl (5i)/-0 2 , Zl (ii), as required by (2.27), 
coincides with q < z ] ) Zl , as defined in (2.33), 

( 2 -43) 

By formula (2.27) for p zi , zi and Lemma 2.2 (i) for p zi — — V ; 2(zi)/V'i( z i) onc infers 
p Zl , Z1 (x) = -p' Zl (x) + q { z l~ zi (x)p Zl {xf + 2iS lPzi (x) 

= - (q(x)p Zl (x) 2 + 2izip Zl (x) - p(x)) + q ( ^ rzi (x)p Zl {xf + 2izip Zl (x) 

= Pzl~zM- ( 2 - 44 ) 
Using (2.29) and (2.22)-(2.24) onc computes 

T(z,p Zl ,q Zl . Zl )f{z,q,p Zl ) 

-(i/2)2r/ 2 -(i/2)W(*(a;i),*(5i)) 1 



x 



•Si-02(zi)-0i(ii) - 2:1 ipi (21)^2(^1) -(^1 7 2:1)^1 (zi )i>i(zi)_ 

(Zi - Zi)lp 2 (zi)V> 2 (^l) -Zi-01 (21)^2(21) + 21^2(^1)^1(^1) 



= r(«, *(«!)). (2.45) 

To check (2.35) one uses Lemmas 2.4 and 2.5, 
H*, *(*!),*(;?!)) 

= f'(2,p Zl ,^ li5l )f (z,q,p Zl ) + t(z,p Zl ,q ZltZl )t'(z,q,p Zl ) 
= [-^(z,p Zl ,q Zuil )U(z,p Zl ,q Zl ) 

+ U(z,p ZuZl ,q ZuZl )t(z,p Zl ,q ZltZl )]t(z,q,p Zl ) 
+ t(z,p Zl ,q ZuZl )[-f(z,q,p Zl )U(z,p, q) + U(z,p Zl ,q Zl )f(z,q,p Zl )\ 
= U(z,p Zl , Zl ,q Zl , Zl )T(z, *(*!),*(«!)) - F(z, *(*i))tf(*,p, g). (2.46) 

Formulas (2.39)-(2.40) follow from (2.32)-(2.34) since by Lemma 2.2 (iii), 

W(*( Zl ), *(2T)) - /Ctt(si)) = -||*(*i)|&. (2.47) 

□ 
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For a general treatment of Backlund (Darboux) and gauge transformations and 
their interrelations we refer, for instance, to [34, Sect. 4.1] and [39]. 
Finally we add a few more facts valid in the NLS_ case. 

Remark 2.7. Assume the NLS_ case p = —q. If T(z) = (v 1 (z), v 2 (z)) T is a 
z-wave function associated with q, then lCT(z,x) is a 'z-wave function associated 
with q and 

T« (z, x) = T(z, x, *(^), /Ctt(zi))T(Z, x) 
= T(z, x, *(zi), /C*(zi))/CT(z, x) 

= KT^(z,x) (2.48) 
is a z-wave function associated with qi\^ (cf. Lemma 2.2(iii)). 

Remark 2.8. Assume the NLS_ case p = —q. 

(i) Take z = Z\ and T(zi) = \&(zi) in Theorem 2.6. Then 

*%>(z u x) = T(z 1 ,x,^(z 1 ),JC^(z 1 ))^(z 1 ,x) =0 (2.49) 
since *(z i ) ± *(zi) = 0. 

(ii) Take z — zi and /CT(zi) = /C*(zi) in Theorem 2.6. Then 

K,¥^{ Zl ,x) = r(zT,a;,*(zi),/C*(zi))/C*(zi,a;) = (2.50) 
since *(zi)- L /C*(zi) = ||\l>(zi)||^ 2 by Lemma 2.2{iii). 

In the NLS_ case, Theorem 2.6 also shows that qi\^ is locally nonsingular when- 
ever q is. More precisely, one has the following result. 

Corollary 2.9. Assume the NLS_ case p = —q and suppose z\ G C. Then, if 
q G I/f oc (M) for some p G [1, oo) U{oo} (respectively, if q G C k (R) for some k G No), 
the NLS- potential qi\^ given by (2.38) also satisfies q^ G if oc (R) (respectively, 
qi? G C*(R)). 

Proof. Since |^i(zi, x)ip2(zi, x) | ||W(zi, a; ) 1 1 < 1/2, one concludes from (2.38) that 

(qil } - 0) e L°°(R). Again by (2.38), q and q^ share the same L\ oc and C k 
properties since by (2.4) (with p = —q) one has d™ij)j(z\, •) G ACi oc (R), j = 1,2, 
whenever d™q G L\ oc (R). □ 

3. J'-Self-Adjointness of NLS_ Dirac-Type Operators 

It is a known fact that the Dirac-type Lax differential expression in the defocusing 
NLS+ case is always in the limit point case at ±oo. Put differently, the maximally 
defined Dirac-type operator corresponding to the defocusing NLS + case (cf. (1.3)) 
is always self-adjoint. Classical references in this context are [32, Sect. 8.6], [45], 
which use some additional conditions (such as real-valuedness and/or regularity) of 
the coefficient q. A simple proof of this fact under most general conditions on q was 
recently communicated to us by Hinton [23] (cf. also [5], [6] and [29] for matrix- 
valued extensions of this result). In this section we show that the analogous result 
holds for Dirac-type differential expressions M(q) in (2.7) in the focusing NLS- 
case, when self-adjointness is replaced by J'-self-adjointness. 
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First, we recall some basic facts about j7-symmetric and j7-self-adjoint operators 
in a complex Hilbert space H (see, e.g., [11, Sect. III. 5] and [21, p. 76]) with scalar 
product denoted by (■, -)n (linear in the first and antilinear in the second place) 
and corresponding norm denoted by || • \\n- Let J' be a conjugation operator in H, 
that is, J' is an antilinear involution satisfying 

(Ju, v) n = (Jv, u) H for all u, v G H, J 2 =I. (3.1) 

In particular, 

( Ju, Jv) H = (v, u) H , u,veH. (3.2) 
A densely defined linear operator S in H is called J'-symmetric if 

S C JS*J (equivalently, if JSJ C S*). (3.3) 

Clearly, (3.3) is equivalent to 

{Ju, Sv) H = (JSu, v) H , u, v e dom(S). (3.4) 

Here S* denotes the adjoint operator of S in H. If S is j7-symmetric, so is its 
closure S. The operator S is called J-seU- adjoint if 

5" = JS*J (equivalently, if JSJ = S*). (3.5) 

Finally, a densely defined, closable operator T is called essentially ^-self-adjoint if 
its closure T is J-self adjoint, that is, if 

T = JT*J. (3.6) 

Next, assuming S to be J'-symmetric, one introduces the following inner product 
(•, •)* on dom( JS*J) = dom(S*J) according to [26] (see also [37]), 

(u, v)„ - ( Ju, Jv) H + (S*Ju, S*Jv) H , u,ve dom(JS*J), (3.7) 

which renders dom(j7 S*J) a Hilbert space. Then the following theorem holds (I 
denotes the identity operator in TL). 

Theorem 3.1 (Race [37]). Let S be a densely defined closed J -symmetric operator. 
Then 

dom(JS*J) = dom(S) ©* ker((S**J) 2 + /), (3.8) 

where ©* means the orthogonal direct sum with respect to the inner product (■, •)*. 
In particular, a densely defined closed J -symmetric operator S is J -self-adjoint if 
and only if 

kcr((S*J) 2 +I) = {0}. (3.9) 

We will apply (3.9) to (maximally defined) Dirac-type operators associated with 
the differential expression M(q) in (2.7) relevant to the focusing NLS_ hierarchy 
and prove the fundamental fact that such Dirac operators are always J'-self-adjoint 
under most general conditions on the coefficient q (cf. Theorem 3.5). 

It will be convenient to make the following NLS_ assumption throughout the 
remainder of this section. 

Hypothesis 3.2. Suppose q £ Ll QC (M) and assume the NLS_ case p = —q. 
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Given Hypothesis 3.2, we now introduce the following maximal and minimal 
Dirac-type operators in L 2 (R) 2 associated with the differential expression M(q), 

D max (q)F = M(q)F, (3.10) 

F 6 dom(Anax(<z)) = {G e L 2 (R) 2 \ G e AC loc (R) 2 : M(q)G e L 2 (R) 2 }, 

D min (q)F = M(q)F, (3.11) 

F e dom(D mi n((?)) = {G G dom(L> max (g)) | supp(G) is compact}. 

It follows by standard techniques (see, e.g., [32, Ch. 8] and [45]) that under Hypoth- 
esis 3.2, D m i n (q) is densely defined and closablc in L 2 (R) 2 and D max (q) is a densely 
defined closed operator in L 2 (R) 2 . Moreover (cf. (2.18)), one infers (see, e.g., [32, 
Lemma 8.6.2] and [45] in the analogous case of symmetric Dirac operators) 



Anin(g) = D max (-q)*, or cquivalently, L> m in(<?)* = D max (-q). (3.12) 

The following result will be the crucial ingredient in the proof of Theorem 3.5, 
the principal result of this section. 

Theorem 3.3. Assume Hypothesis 3.2. Let N(q) be the following (formally self- 
adjoint) differential expression 

N(q)=i(£ 'A (3.13) 
V 1 dx/ 

and denote by D max (q) the maximally defined Dirac-type operator in L 2 (R) 2 asso- 
ciated with N(q), 

D ma MF = N(q)F, (3.14) 
F e dom(A nax (g)) = {Ge L 2 (R) 2 | G e ACi oc (R) 2 ; N(q)G e L 2 (R) 2 }. 

Then one infers: 

(i) The following identity holds 

M(-q)M(q) = N(q) 2 . (3.15) 

(ii) Let U q = U q (x) satisfy the initial value problem 

U' q =(^- tyu q , U q (0) = I 2 . (3.16) 

Then {U q (x)} x& R is a family of unitary matrices in C 2 with entries in ACi oc (R) fl 
L°°(R) satisfying 

U- 1 N(q)U q = i-^L 2 . (3.17) 

(in) LetU q denote the multiplication operator with U q () on L 2 (R) 2 . Then D max (q) 
is unitarily equivalent to the maximally defined operator in L 2 (R) 2 associated with 
the differential expression i-^.li, 



Uq'D^iqpq = ( l — h , (3.18) 

aX / max 

= H h2 (R) 2 = {F e L 2 (R) 2 | F G ACi oc (R) 2 ; F' £ L 2 (R) 2 }. 



domf ( i—I' 
\ \ dx 
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Moreover, 

U- 1 D mm (-q)D max (qp q = (-^2) (3.1!)) 
(loin ( ( -^I 2 ) )=H 2 - 2 (R) 2 



<( 



dx 2 2 



{F e L 2 (R) 2 \F,F' G ACi oc (K) 2 ; G £ 2 (R) 2 }. 



Proof. That 7V(<?) is formally self-adjoint and M(-q)M{q) = N(q) 2 as stated in (i) 
is an elementary matrix calculation. 

To prove (ii), we note that the initial value problem (3.16) is well-posed in the sense 
of Caratheodory since q e L 1 1 oc (E) (cf., e.g., [22, p. 45-46]) with a solution matrix 
U q with entries in ACi oc (M.). Moreover, for each x G M, U q (x) is a unitary matrix 

in C 2 , since U' q = —B(q)U q , with B(q) = skew-adjoint. Thus, the entries 

U q ,j,k, 1 < j, k < 2 of U q (as well as those of Uq 1 ) actually satisfy 

U qJ , k G ACi oc (R) n L°°(M.), 1 < j, fc < 2. (3.20) 

Next, fix F G AC loc (M) 2 , such that Z^F G i/ 1,2 (R) 2 . Then 

= ? ;^F + zC/ 9 (c/- 1 B( (Z r)F 

= (3.21) 

where we used the fact that {Uq 1 )' = U q 7 1 B(q)* . Thus, (ii) follows. 

Moreover, by (3.21) one concludes dom(5 max (g)) = UqH 1 ' 2 ^) 2 by (3.20) and the 

fact that U q is unitary in C 2 . This proves (3.18). 

Clearly (i) and (ii) yield the relation 

U^M{-q)M{q)V q = -^h- (3.22) 
Thus, (3.19) will follow once we prove the following facts: 

(i) U q F G L 2 (R) 2 if and only if F G i 2 (R) 2 , 

(ii) J7 9 F e ACioc(M) 2 if and only if F G ACi oc (M) 2 , 

(iii) M(q)U q F £ L 2 (R) 2 if and only if F' G L 2 (R) 2 , 

(iv) M(q)U q F G ACi oc (M) 2 if and only if F' G ACioc(K) 2 , 

(v) M(-q)M(q)U q F G P 2 (R) 2 if and only if F" G L 2 (R) 2 . 

Clearly (3.23) and (3.27) hold since U q is unitary in C 2 . (3.24) is valid since 
Uq.j^u, U~j k G ACioc (R) n L°° (R) , j, k = 1,2. Next, an explicit computation yields 

M(g)tf,F = i( U T f' lf ^ U rf' 2f i), F=(f 1 ,f 2 ) T . (3.28) 

\ — U q,2,lJl - Uq,1,2)l) 



(3. 


.23) 


(3. 


.24) 


(3. 


.25) 


(3. 


.26) 


(3. 


.27) 



Introducing 
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one infers Vq,j,k,V~J k 6 ACi oc (R) n L°°(K), j, k = 1, 2 and 

F- 1 Af( <Z );7 (? F = z(/{,-.^) T , (3.30) 
and hence (3.25) and (3.26) hold. This proves (3.19). □ 

Remark 3.4. We note that by (3.18), D max (q), the maximally defined, self-adjoint 
operator in L 2 (R) 2 associated with the 2x2 matrix-valued differential expression 
N(q) has a purely absolutely continuous spectrum that equals R, 

<r{D max {q)) = MAnaxfa)) = R - ( 3 ' 31 ) 

( We refer to Section 5 for a discussion of various spectral decompositions. In the 
present context we just note that a(T) and cr ac (T) denote the spectrum and ab- 
solutely continuous spectrum of a self-adjoint operator T in a separable complex 
Hilbert space TL.) 

The principal result of this section then reads as follows. 

Theorem 3.5. Assume Hypothesis 3.2. Then the minimally defined Dirac-type 
operator D m i n (q) associated with the Lax differential expression 

M(q) = i(& ~f) (3.32) 

introduced in (3.11) is essentially J -self-adjoint in L 2 (M.) 2 , that is, 

Amn(<z) = JD min {q)*J, (3.33) 
where J is the conjugation defined in (2.9). Moreover, 

Anin(g) - Anax(g) (3-34) 

and hence D max (q) is J -self-adjoint. 
Proof. We first recall (cf. (3.12)) 

Anin(<7)* - Anax(-<Z) (3.35) 

and also note 

JD mm (-q)J = D mm (q). (3.36) 

Since -D m in(g) is closed and J'-symmetric (cf. (2.19)), its J'-self-adjointness is equiv- 
alent to showing that (cf. (3.9)) 

kcr (D min (q)*JD min ( q yj + I)= kcr(A„ax(-g)Anax(g) + /) = {0}. (3.37) 

Since by (3.19) D max (— q)D max (q) is unitarily equivalent to {—d 2 /dx 2 l2) m a,x > 0, 
one concludes D max (—q)D max (q) > and hence (3.37) obviously holds. The fact 
(3.34) then follows from (3.33) and (3.35) since 

Aninfa) = JD min (q)*J - JD max (-q)J = D max (q). (3.38) 

□ 

As mentioned in the introductory paragraph to this section, Theorem 3.5 in the 
^-self-adjoint context can be viewed as an analog of the result of the corresponding 
(self-adjoint) Dirac operator relevant in the defocusing NLS+ case of the nonlinear 
Schrodingcr equation. 
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4. Constructing L 2 (R) 2 -Wave Functions for J'-Self- Adjoint 
Dirac-Type Operators 

In this section we discuss how to construct L 2 (R) 2 -wave functions for non-self- 
adjoint (but J'-self-adjoint) Dirac-type operators associated with the Lax differen- 
tial expression for the NLS_ system. 

By Remark 2.8, in order to obtain a nonzero zi-wave function associated with 
qii , we have to apply the transformation matrix T(zi, ^>(zi), K$>(zi)) to a zi-wave 
function 3>(zi) associated with q that is linearly independent with the original Z\- 
wave function \I/(zi) associated with q. Similarly, in order to obtain a nonzero 
zT- wave function associated with qi\\ we have to apply the transformation matrix 
r(zl, ^(zi), K.^(zi)) to a zT- wave function /C$(zi) associated with q that is linearly 
independent with the original zT-wave function /C*(zi) associated with q. The 
function $(zi) is constructed as follows. 

Let *(z) = (V>i(z), ip2{z)) J , z £ C, be a z-wave function associated with goiil 
and introduce 

n z = {x 6R|^(«,jf)^(«,i) ^ 0}. (4.1) 

Next, consider 

**{z,x) = {l/2){i, 2 {z,x)-\-^ 1 {z,x)- 1 ) T , xen z , (4.2) 

such that 

W{**(z,x),*{z,x)) = 1. (4.3) 
Let xo, x £ Q z such that [xo, x] C Q, z and define 

R(z,x,x ) = - 1 - [ X dx '( f X '\ + .f^A (4.4) 

$(z,x) = * # (z,x) + i?(z,x,x )*(z,x), [x ,x]CQ z . (4.5) 

Using (2.4), we have (*#)' = U(p,q)V# - R'V. Thus, W($, ) = 1 and $(z) is a 
z-wave function associated with q, since 

$' = + R'^J + Rx&> = U^* - R'V + R'V + RU^> = {/$. (4.6) 

We note that f- 1 ** = 1 implies /C*(z)*(z) ± $(z) = /C*(z). 

Thus, if F(z, *(zi),/C*(zi)) is defined using <J>(zi) as in (2.40), then the zi-wave 
function $^(zi) associated with as prescribed in Theorem 2.6, is computed 
as follows 

( Zl> x ) = r(zi, x, *(zi), /C*(zi))$(z l5 x) 

= Im(z 1 )||*(z 1 ,x)||- 2 /C*(z 1 , a; ). (4.7) 

Moreover, by Remark 2.7, /C$i^(zi) is computed as 

(zx, x) = F(zl, x, /C*(z 1 ))/C$(z 1 , x) 

= -Im(z 1 )||*(z 1 ,x)||- 2 2 *(z 1 ,x). (4.8) 
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By (4.7) and (2.40), for each z E C, the z-wave function &i\\z) associated with 
qi\^ (constructed using the z-wavc function $(z) associated with q) is computed by 

*W(z, x) = T(z, x, *( Zl ),lC*(zi)Mz, x) 

= -{i/2)(z - Z!)$(z,x) + aO-'-fcCe, a;) 

= -(i/2)(z - zi)$(z, x) - W(*(zi,x),*(z, x))$W ( ZliX y (4.9) 
Formulas (2.38) and (4.7) now imply 

q^>{x)=q{x)+A4^l i {z 1 ,x)^ 1 {z u x), (4.10) 
where *$>{z,x) = (<j>^ Zl (z, x), ^ (z, x)) T , V(z,x) = V> 2 (*, x)) T . 

Remark 4.1. We emphasize that while R(z,x,x ) in (4.4), and ftence $(z,x) m 
(4.5), in general, will have singularities on R, the formulas (4.7)-(4.10) are well- 
defined for all x E R. 

The next hypothesis will be crucial in our attempt to construct z\- and zi-wave 
functions in L 2 (R) 2 associated with the Dirac-type differential expression M{q^\ 

Hypothesis 4.2. Suppose q E L 1 1 oc (IR), assume the NLS_ case p = —q, and let 

z E C. Suppose ty(zo) to be a z -wave function associated with q that satisfies the 
condition \\^(z 0l -)\\^2 E L 2 (R), that is, 

/oo 
dx\\*(z ,x)\\ci < oo. (4.11) 

If a z -wave function ^(zo) associated with q satisfies condition (4.11), we will 
henceforth say that &(zo) satisfies Hypothesis 4.2 at z . 

Remark 4.3. Assume Hypothesis 3.2 and suppose that ^(z) satisfies Hypothe- 
sis 4-2 at z. Then, 

(i) by Holder's inequality, V(z) g L 2 ((-oo, R\f U L 2 ([R, oo)) 2 for all ReR. 

(ii) K$>(z) satisfies Hypothesis 4-2 at~z by Lemma 2.2 (Hi). 

Remark 4.4. Assume Hypothesis 3.2 and let AeR. Then all X-wave functions 
^(A) associated with q satisfy ||^(A, a;)||c 2 = C W {independently of x E R). In 
particular, ||*(A, -)||c 2 i ll*(A, • ) 1 1 c ^ ^ ^ 2 ([^i R EM., and hence there exists 

no \-wave function associated with q that satisfies Hypothesis 4-2 at A £l. 

The principal result of this section then reads as follows. 

Theorem 4.5. Assume Hypothesis 3.2. Let z\ E C\R and suppose that the z\-wave 
function ^(zi) associated with q satisfies Hypothesis 4-2 at Z\. Let q^ be given 
by (2.38). Then Z\ and ~z~\ are eigenvalues of the maximal Dirac-type operator 
fmax(?z!') associated with M(gi^) of geometric multiplicity equal to one. The 
corresponding eigenf unctions <f>i^(zi) and /C$i^(zi) are given by (4.7) and (4.8), 



respectively, that is, one has 

^(z 1 ),K^{z 1 ) E dom(^ max (gW)), (4.12) 

D msK (q^)^(z 1 ) = z 1 ^(z 1 ), (4.13) 

O™^')^^) = zllC*>£>(z 1 ). (4.14) 
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Proof. Indeed, using (4.11) at z\ one obtains 

/oo 
dx\\*{z u x)\\^\\*{z u x)f c , <<x>. (4.15) 
-oo 

In order to show that z\ has geometric multiplicity equal to one as an eigenvalue of 
D max (q z ^), we next assume that ^(z lr ) G L 2 (R) 2 and ^(z lr ) G L 2 (R) 2 
are linearly independent ^i-wave functions associated with qi\\ Then clearly 
^(z u -),^(z u -) G dom(D ma 4q { z i ) )) and 

W{-fi\Hz lr ),¥ z \Hz lr ))eL 1 (R). (4.16) 

However, since by (2.5), W{§z) (zi,x),&z) {z\,x)) is constant with respect to x G 
R, (4.16) represents a contradiction and hence z\ has geometric multiplicity equal 
to one. The analogous arguments apply to zj. □ 

The argument that -D m ax(<zii^) has only eigenvalues of geometric multiplicity 
equal to one applies of course in complete generality to any D{q) with q G L 1 1 oc (R). 

Next we show that condition (4.11) is preserved under iterations, a fact of great 
relevance in connection with the multi-soliton solutions relative to arbitrary back- 
grounds discussed at the end of Section 6. 

Lemma 4.6. Assume Hypothesis 3.2 and let z\,z 2 G C. Fix a zi-wave function 
^{zi) and a z 2 -wave junction $(22) associated with q. Using ty(zi), construct 
the NLS_ potential qi]^ by formula (2.38) and consider the transformation matrix 
T(z2, *(2;i),/C*(zi)) given by formula (2.40) for z = z 2 . Let 3> z \\z 2 ) be the z 2 - 
wave function associated with q z \^ as in (4.9). If ®(z 2 ) satisfies Hypothesis 4-2 at 
z 2 , then ®iWz 2 ) satisfies Hypothesis 4-2 at z 2 . 

Proof. Without loss of generality we may assume Im(0i)Im(z2) > 0. Formulas 
(2.12), (4.7), and (4.9) imply 

= \\Z2- z 1 \ 2 \\Hz 2 )\\l 2 + \\zl- z 1 | 2 ||*(z 1 )||c 2 4 ||/C*(z 1 ))vI/(z 1 ) ± <I'(z 2 )||^ 

-iRc(((z 2 -z 1 )<I>(z 2 ),(zr-z 1 )||*(z 1 )||- 2 2 /CvI/(z 1 )vI/(z 1 ) ± $(z 2 )) c2 ) 
= \\z 2 - zi| 2 ||*(^)||^ + \\z- 1 -z 1 \ 2 \\^(z 1 )\\^(z 1 )^(z 2 )\' 2 



1 

2 

1 



+ -\\*(z 1 )\\^(z 1 )^(z 2 )\ 2 Rc((z 2 - Zl ){zi - Zl )) 



- A \z 2 zi| 2 ||$(z 2 )|| 2 2 + ii*^)!!- 2 !*^)^^^)! 2 ^^)^^) 
1 



> i |z 2 ~z 1 | 2 ||$(z 2 )|| 2 2 . (4.17) 

Thus, if j x oc dx\\^(z 2 ,x)\\- 2 < 00, then f™ x dx \\^{z 2 , x) < 00. □ 

It will be shown in Remark 6.1 that for any z G p{D Yaa-yi ( y q^ > )) all but two z-wave 
functions of £> m ax(?zi^) satisfy Hypothesis 4.2 at z. 
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5. Some Spectral Properties and the Existence of 
Weyl-Titchmarsh-type Solutions for J'-Self-Adjoint 
Dirac-Type Operators 

The principal purpose of this section is to establish the existence of Weyl- 
Titchmarsh-type solutions for formally J'-self-adjoint Dirac differetial expressions 
M(q) associated with the focusing NLS_ case. The latter are well-known to ex- 
ist in the case of self-adjoint Dirac operators (in particular, they are well-known 
to exist in the context of the defocusing NLS+ equation) and are known to be a 
fundamental ingredient in the spectral analysis in the self-adjoint context (cf., for 
instance, [31, Chs. 3, 4,]). As far as we know, no such result appears to be known 
in the general ^7-self-adjoint case studied in this section. Along the way we also 
collect some results concerning spectral properties of D max (q). 

Thus, assuming the NLS_ case and hence the basic Hypothesis 3.2 throughout 
this section, we adopt the simplified notation (cf. Section 3) 



D(q) = D min (q) = D max (q), q G L\ oc (R) . (5.1) 

We also denote by I 2 the identity operator in L 2 (R) 2 (as well as in C 2 ). Moreover, 
we find it convenient to introduce the following notations, 

L? oc ([-oo,oo)) = {/: 1R — > C is measurable | / G L 2 ((-oo,R]) for all R G K}, 

Lf oc ((-oo, oo]) = {/: R-> C is measurable | / G L 2 ([R,oo)) for all R G M.}. 

(5.2) 

We start with the following auxiliary result (the variation of parameters formula) . 

Lemma 5.1. Assume Hypothesis 3.2 and let (z, Xq) £ CxM. Let ^i(z) and ^>2{z) 
be linearly independent z-wave functions for M(q) defined on [xq, 00) and denote by 
E(z, x) = [*i(-2, x), ^2{z, x)] the 2x2 fundamental matrix solution of M(q)E(z) = 
zE{z). Assume the Wronskian of ^1 and ^2 satisfies W(^i(z, x), ^2(2, x)) = 
dct(S(z,x)) = 1 for some (and hence for all) x G [xo, 00). Moreover, suppose B G 
L 1 1 oc ([x ,oo)) 2x2 . Then <f>(z,-) G AQ oc (R) 2 satisfies M(q)$(z) = z$(z) + 5$(z) if 
and only if 

®(z,x) = E(z,x)C + I dx'E(z,x)E(z,x')~ 1 B(x')^(z,x'), x > x (5.3) 

J x 

for some C = (ci, c 2 ) T G C 2 independent of x. Moreover, &(z, x) = for all x > x 
if and only if C = 0. 

Proof. The computation 

M(g)($(x) - r ( ix'S(x)S(x')" 1 S(x')*(x')) 

J Xo 

= z<f>(x) + B(x)$(x) - M(q) (z(x) J dx' S(x') _1 B(x')3>(x')^ 

= z ($(x) - £ dx' E(x)E(x')- 1 B(x')<i>(x')^ (5.4) 

shows that $ satisfies M(q)$(z) = z$(z)+B$(z) if and only if $ satisfies (5.3) since 
S is a fundamental matrix for the first-order linear differential system M(q) i S(z) — 
zty(z). That <j>(z,x) = for all x > x if and only if C = follows upon iterating 
the Volterra-type integral equation (5.3) in a standard manner. □ 
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Next, we find it convenient to recall a number of basic definitions and well-known 
facts in connection with the spectral theory of non-self-adjoint operators (we refer 
to [11, Chs. I, III, IX], [21, Sects. 1, 21-23], and [38, p. 178-179] for more details). 
Let S be a densely defined closed operator in a complex Hilbert space TL. Denote 
by B(TL) the Banach space of all bounded linear operators on TL. The spectrum, 
<j(S), point spectrum (the set of eigenvalues), & P (S), continuous spectrum, a c (S), 
residual spectrum, a T (S), approximate point spectrum, cr a (S), essential spectrum, 
a e (S), field of regularity, tt(S), resolvent set, p(S), and A(S) are defined by 



<7(S) = {\eC\(S-XI)- 1 eB(H)}, (5.5) 

<j p {S) = {X G C| kcr(S- XI) ± {0}}, (5.6) 
a c (S) = {Xe a(S) | ker(S - XI) = {0} and ran(S - XI) is dense in TL 

but not equal to TL}, (5-7) 

a T (S) = {X G C | kcr(S' - XI) = {0} and ran(5 - XI) is not dense in TL}, (5.8) 
&a,(S) = {X G C | there exists a sequence {/„}„eN with ||/ n ||« = 1, n G N, 

and lim - XI)f n \\ H = 0}, (5.9) 

n^oo 

o~ e (S) = {A £ C | there exists a sequence {/ n } n eN s -t. {/ n }nGN contains no 

convergent subsequence, ||/ n ||w — 1, n G N, and lim ||(5 — \I)f n \\u = 0}, 

n— *oo 

(5.10) 

tt(5) = {z G C I there exists k z > s.t. \\(S - zl)u\\ n > k z \\u\\ H 

for all u G dom(5)}, (5.11) 

p(S) = C\a(S), (5.12) 

A(S) = {z G C | dim(kcr(5' - z)) < oo and ran(S' - z) is closed}, (5.13) 
respectively. One then has 

a(S) = (T P (S) Ua c (S) U a T (S) (disjoint union) (5.14) 

= <7 P {S) Ua e (S) Ua T (S), (5.15) 

<J C (S) C <j e (S)\(a p (S) U a T (S)), (5.16) 

a T (S) - a p (S*)*\a p (S), (5.17) 
<7*(S) = {A G C | for all e > 0, there exists 7^ f s 

s.t. ||(5-A)/ £ || w <s||/ £ || H } (5.18) 

= C\ir(S), (5.19) 

t7(5)\<7 B (5) C o- r (5), a(S)\a a (S) is open, (5.20) 
cr e (5) = {A G C I there exists a sequence {g n }nen s -t- w-linin^oo g n = 0, 

\\9n\\n = 1, n € N, and lim ||(5 - A/) 5 „|| w - 0} (5.21) 

-co 

= C\A(S), (5.22) 

Ce(-S) Q &a(S) C ct(S') (all three sets are closed), (5.23) 

C ^(S 1 ) C A(5) (all three sets are open). (5.24) 

Here uu* in the context of (5.17) denotes the complex conjugate of the set u C C, 
that is, 

uj* = (X G C| A G u}. (5.25) 
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a(S) 


= <t p (S)Uct c (S) 




= a p (S)Ua e (S), 


<Jr(S) 


= 0, 




= a p (S*y, 


a a (S) 


= a(S), 


tt(S) 


= P(S), 



For future reference we note that the sequence {f n }neN hi (5.10) (and the sequence 
{<?n}neN in (5.21)) is called a singular (or Weyl) sequence of S corresponding to 
A. We also note that there are several versions of the concept of the essential 
spectrum in the non-self-adjoint context (cf. [11, Ch. IX]) but we will only use the 
one in (5.10) (respectively, (5.21)) in this paper. 

In the special case where S is J'-self-adjoint one obtains the following simplifi- 
cations (cf. [11, p. 118], [21, p. 76]): 

(5.26) 
(5.27) 
(5.28) 
(5.29) 
(5.30) 
(5.31) 

whenever 5* is j7-self-adjoint. Note that 7r(S) = may occur for ^-symmetric op- 
erators (see [37]) in sharp contrast to the case of densely defined, closed, symmetric 
operators T, where any nonreal number is in the field of regularity n(T). 

Returning to the NLS_ case, we next recall an elementary but useful consequence 
of (2.16) and (2.18). 

Lemma 5.2. Assume Hypothesis 3.2. Then 

JCD(q)K, = -D(q), (5.32) 
<T3D(q)a 3 = D(q)* = D(-q). (5.33) 

Consequently, 

a(D(q)r=a(D(q))=a(D(- q )). (5.34) 

Moreover, 

if z a e a p (D(q)) and D(q)F = z a F for some F £ dom(D(g)), 

then D(q)K,F = z^JCF. (5.35) 

Proof. Relations (5.32) and (5.33) are clear from (2.16), (2.18), (3.10), (3.12), and 
(5.1). Since (D(q)* - zh)- 1 = [(D{q) -zh)- 1 }* , a(D(q)*) = a(D(q))* and hence 
(5.34) follows from (5.33). Finally (5.35) is clear from (5.32) and K 2 = -I 2 - □ 

Next we introduce the basic hypothesis to be assumed for the remainder of this 
section. 

Hypothesis 5.3. Suppose q G L 1 1 oc (M), assume theNLS- case p = —q, and suppose 
that the (J -self-adjoint) operator D(q) has nonempty resolvent set, p(D(q)) ^ 0. 

Recalling the standard notation 

nul(T) = dim(kcr(T)), (5.36) 

def(T) = dim(ran(T)- L ) = dim(ker(T*)) = nul(T*), (5.37) 

where T denotes a densely defined closed operator in Ti., we can state the following 
fundamental result, establishing the existence of Weyl-Titchmarsh-type solutions 
for J'-self-adjoint Dirac-type operators relevant to the NLS_ case. 
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Theorem 5.4. Assume Hypothesis 5.3 and pick z G p{D(q)). Then there exist two 
unique (up to constant multiples) linearly independent z-wave functions ^-(z, •) 
and ^+{z, ■) associated with q satisfying 

*_(*,•) €£^([-00,00)), *+(*,•) e£L((-°°.°°]). ( 5 - 38 ) 

\\9-(z,-)\\£ eLf oc ((-oo,oo]), ||vD, + (^, e if oc ([-oo, oo)), (5.39) 
lim ||* ± (*,aO|| C 2= lim \\*±(z,x)\\^ = 0, (5.40) 



sup 



| r ds||*_(*,aOI&) QHds||* + (z,aO|& 



< oo. (5.41) 



The corresponding z-wave functions K,^±(z,x) associated with q satisfy (5.38)- 
(5.41) with z replaced by z. 

Proof. We prove the existence of the two z-wave functions following the lines of 
[11, Sect 10.4]. To this end we introduce the operators 

D mgtK (q;-^)F = M(q)F, 

F e dom( J D max ( g; -co)) = {Ge £ 2 ((-oo, 0]) | G e AGi oc ((-oo, 0]), (5.42) 

M(q)G G L 2 ((— oo, 0])}, 

D mia (q; -^)F = M(q)F, 

F G dom(D min (q; -oo)) = {G G dom(D max (q; -oo)) | G(0) = 0; (5.43) 

supp(G) C (—00,0] is compact}, 

Anax(?; oo)F = M(q)F, 

F G dom(£» max (g; oo)) = {G G L 2 ([0, oo)) | G G AGi oc ([0, oo)), (5.44) 

M(q)G Gi 2 ([0,oo))}, 

£) min ( g ;oo)F = M(g)F, 

dom(U min ( g ; oo)) = {G£ dom(U max (g; oo)) | G(0) = 0; (5.45) 
supp(G) G [0, oo) is compact}. 
In close analogy to [11, Theorem III. 10. 20] one can prove that for all z G p(D(q)), 



dcf(Anin(<z) - zh) = def(Anin(<Z; -oo) - zh) + def( J D min ( g; oo) - zl 2 ) - 2. (5.46) 



Since D(q) = D m \ n (q) is J'-sclf-adjoint and z G p(D(q)), we necessarily have (see, 
e.g., [11, Theorem III.5.5]) that 



dcf(D min (g) - zl 2 ) = Aci{D(q) - zl 2 ) = 0. (5.47) 

Thus, 

def(D min (g; -oo) - zl 2 ) + dei(D min (q; oo) - zl 2 ) = 2. (5.48) 
We claim that the only possibility is 



def(Z? min (g; -oo) - zl 2 ) = dcf(L> min (g; oo) - zl 2 ) = 1. (5.49) 
Indeed, arguing by contradiction, we assume, for instance, that 



def( J D min (g; -oo) - z a I 2 ) = 0. (5.50) 

This implies 

nul(D max (-g; oo) - zqI 2 ) = 2 (5.51) 
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since (D m i n (q); oo)* = D max (— q; oo). Thus all zo-wave functions for the NLS_ 
potential — q are in Lf oc ((— oo, oo]) 2 . According to Remark 4.4, this is clearly im- 
possible for zq G K fl p(D(—q)). Next we show that this is impossible also for 
zo G (C\M) n p(D(—q)). Using Lemma 2.2 (hi), we can simplify notations and re- 
place zq by zo without loss of generality. Moreover, to avoid confusion with the 
change q — ► —q and the corresponding change for q^J , we simply use q instead of 
—q in the proof of (5.38) below. 

To this end, we hx z G p(D(qj) with Im(z ) ^ and *i(^o), ^2(^0) two lin- 
early independent Zo-wave functions associated with the background potential q. 
The latter are in L 2 ([0, oo)) 2 by hypothesis (5.51). Then for any zo-wave function 
vp(zo) = (i>i(zo), ip2{zo)) T associated with the NLS- potential q, one infers (cf. 
(4.7) and (4.10)) that $<i \zo) = lm(z ))C^ (z )\\^ {z )\\ C 2 is a zo-wave function 
associated with the NLS_ potential 



q + AIm(z )ipi(zo)Mzo)\\^{zo)\\ C 2. 
Thus $zJ(zo) satisfies 

z $W(z ) = M(q^)^(z ) = M(g)*£>(*o) - B(z )^(z ), 

where 





B(z , x) = 4ilm(z )\\y{z ,x)\y 
belongs to L a 



ipi(z ,x)ip 2 (zo,x) 
tpi(z ,x)ip2(zo,x) 



(5.52) 
(5.53) 
(5.54) 



2x2 . In particular, 

esssup a . eK ||B(2;o,a;)||c2x2 < 2|Im(z )|- 



(5.55) 

Since no confusion can arise we occasionally suppress the explicit zo-dependence in 
the calculations below. The variation of parameters formula (5.3) then yields the 
following for the fundamental system of solutions 2 = [^1,^2], W(^i,^2) = 1, 
associated with q and zo: 

= Z(x)C+ f dx'S(x)S(x')- 1 i?(x')<i>i 1 ) (x') (5.56) 

J a 

= c 1 * 1 (x)+c 2 * 2 (x)+S(x) f dx'Zix'^Bix'^ix'), x>a>0, 

J a 

where C = (ci, c 2 ) T 7^ depends on a and we assume W(^i, ^2) = 1 according to 
Lemma 5.1. Hence one obtains 

||S(x)|| C 2x2 <2 1 / 2 max(||* 1 (x)|| C 2,||* 2 (x)|| C 2). 
Moreover, writing fyj = (ipij, ip2,j) , j — 1,2, one infers 

IIV'fcj(2o,-)llL 2 ([a,oo)) < K(a), j,k= 1,2 
for some constant K{a) > with 



lim K(a) = 0. 

a'loo 



(5.57) 
(5.58) 
(5.59) 



Hence, 



rdx , iis(x , )- i iic2x2 < 2 r dx' 

J a J a 



max 



v>2,2(z') 

-MX') 



Mx') 



< 4K(a) 



C 2 / 

(5.60) 
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Thus, (5.55)-(5.60) yield for x > a, 
\\*$(x)\\<? < |ci|||*i(x)|| c2 + \c 2 \\\^ 2 (x)\\ c , (5.61) 

(rx \ 1/2 

J dx'\\$W(x')\\ 2 c2 ) . 

Squaring (5.61) and integrating the result from a to x, one estimates 

[ X dx'\\$W(x')\\ 2 c2 

<3| C1 | 2 r^'n^^oii^+s^i 2 r dx'w* 2 {x')wi 2 

J a J a 

+ 96|Im(z )| 2 ^(a) 2 f dx' f dx" max (||*i(a/)|| c2 , ||*2(a/)||<?) || 

J a J a 

<6K(a) 2 (\ Cl \ 2 + \c 2 \ 2 ) 
+ 96\lm(z )\ 2 K(a) 2 f dx" f dx' max (||*i(a;')||c>, Il*2(a/)||c0 KVOfc 

Ja Jx" 

<6K(a) 2 (\ Cl \ 2 + \c 2 \ 2 ) 
+ 96\Im(z )\ 2 K(a) 2 f dx" \\^ {x")\^ 2 f dx' max (||*i(a/)||<?, Il*2(z')llc0 

J a J a 

< 6K(a) 2 (\ Cl \ 2 + \c 2 \ 2 ) +192|Im(z )| 2 if(a) 4 f dx" \\^{x")f c2 . (5.62) 

J a 

Here we applied the Fubini-Tonelli theorem to the integrand 

ma X (||vI/ 1 ( a; 0|| c2 J|* 2 ( a; ')llcOII < f > i 1 o ) (^)||J 2 X[a,,'](2: , ')>0 (5.63) 
(xa the characteristic function of the set A C t) to prove equality of the iterated 
integrals f x dx' J x dx" ■ ■ ■ and J x dx" J x „ dx' • • • in (5.62). Hence, if one chooses 
a large enough (such that 192|Im(z )| 2 ^(a) 4 < 1), then 

[1 - 192|Im(z )| 2 if(a) 4 ] f dx' \\$<£ (z ,x')\\l 2 < 6K(a) 2 (\ Cl \ 2 + \c 2 \ 2 ), (5.64) 

J a 

and thus, $^ } (z ,-) G L 2 ([a, oo)) 2 . Since H^i^o, -)|| C 2 = H*( z o, Ollc* this con - 
tradicts the assumption ^S(zq, •) G L 2 ([0, oo)) 2 . This proves (5.49). 

Finally, if $-(z) and *+(z) satisfying (5.38) were linearly dependent then it 
would follow that z G a p (D(q)) by (5.38), contradicting the initial assumption 
z G p(D(qj). Summing up, (5.49) implies existence and uniqueness (up to constant 
multiples) of ^f±(z) satisfying (5.38). 

To prove (5.39) we assume without loss of generality that 

W(V-(z,x),V + (z,x)) = l, iGl. (5.65) 

Then one computes 

1 + |(*_(*, x), *+(*, x))c I 2 = ||*+(*, x)|| 2 2 ||*_(z, x)\\ 2 c2 > 1 (5.66) 
and hence, 

\\**{z,x)\\£<\\* ± {z,x)\\c>, ^GK. (5.67) 

Thus, H*.^,.)!^ 1 G £L((-oo,oo]). The fact that ||tt + (z, -)\\£ G L? oc ([-oo, 00)) 
in (5.39) is proved analogously. 
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By (5.38) and (5.39) one infers 

liminf ||*±(z,a;)|| C 2 =liminf ||* T (z,a;)||- 2 1 = 0. (5.68) 

x— »±oo £^±00 

To prove (5.40) one first integrates (2.17) to obtain 

\\* ± (z,x 2 )f c2 - \\*±(z, Xl )\\ 2 c2 = 2Im(z) r dx [|Vi,±M| 2 - hfe,±(*,aO| 2 ]- 

(5.69) 

Thus, lim^^ioo ||\I>±(2, a;)||c 2 exist and hence equal by (5.68). By (5.67) one then 
infers 

lim \\^ T (z,x)\\^=0. (5.70) 

What remains to be proved is (5.41). To this end consider the Green's function 
for the Dirac-type operator D(q) . It can be written in terms of the z-wave functions 



^_{z,x) = (i/j 1 .-(z,x),ip 2 ,-{z,x)) T , V+(z,x) = {ip lt+ {z,x),^2,+ {z,x)) T , 

(5.71) 

whose existence is guaranteed by (5.38) and whose Wronskian is still assumed to 
satisfy (5.65). More precisely, define the 2x2 matrix-valued function on ffi. 2 



*+(*,&)*-(*, !/) T , x < y, 
V-{z,x)V + (z,y) T , x>y, 



G{z,x,y)=-{^r: ! Zrr: y w _z y : ^ 



Explicitly, 



G(z,x,y) = - < 



ipi,+ (z, x)ipi-(z, y) ipi,+ (z, x)ip2,- {z, y) 

K ip2,+ (z, x)ip!-(z, y) ip2,+ {z, x)ip 2 ,- (z, y) / 

i>i- {z, x)tp li+ (z, y) i>i,-(z, x)ip 2 ,+ {z, y) 

ip 2 - (z, x)i/>i l+ (z, y) ip2-(z, x)ip 2 ,+ (z, y) 



x<y, 

(5.73) 

x>y. 



To prove that (5.72) indeed represents the Green's function associated with the 
j7-self-adjoint operator D{q) one can argue as follows. One introduces a densely 
defined operator R(z) in L 2 (R) 2 by 

(R(z)F)(x) = [ dx'G(z,x,x')F(x'), z G p(D(q)), x G R, (5.74) 

Js. 

F G dom(i?(z)) = {G G i 2 (R)) 2 | supp(G) is compact}. (5.75) 

By inspection one infers that 

R(z)F G AC loc (R) 2 n L 2 (R) 2 and M(q)R(z)F G L 2 (R) 2 . (5.76) 

Thus, R(z) maps L 2 (R) 2 -elements of compact support into the domain of D{q). 
Moreover, an explicit computation shows that 

(M(q) - z)[R(z)F - (D(q) - zy'F] 

= {D{q) - z) [R(z)F - (D(q) - z)- 1 F] = 0, F G dom(i?(z)). (5.77) 

Since by hypothesis, z 6 p{D(q)), (5.77) implies 

R{z)F = (D(q) - z)~ 1 F, F G dom(i?(z)). (5.78) 
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Hence R{z) extends boundedly to all of L 2 (R) 2 and its closure coincides with the 
resolvent (D(q) — z)^ 1 of D(q). Thus, z G p(D(q)) is equivalent to the boundedness 
of the operator 

L 2 (K) 2 3 F(x) f G{z, x, y)F(y) dy (5.79) 

JR. 

= *+(*,*) [ y-(z,y) T F(y)dy 

J — OO 

pOO 

+ *_(*,&) / * + (z,y) T F(y)dy 

J X 

in L 2 (R). Taking z G /£>(£>(«)) and F = (/, 0) T and F = (0, /) T it follows that the 
operators 



3 f(x) >-^tp j!+ (z,x) 



I ip k -(z,y)f(y)dy 

J —oo 



fOC 

+ iP jt _(z,x) ij kt+ (z,y)f(y)dy, j,k = 1,2 (5.80) 

J X 

are bounded in L 2 (R). The last statement implies the relations (in fact, it is 
equivalent to them, cf. [35] and Lemma 6.2) 



sup 



a: 



\ip k ,-(x)\ 2 dx 



/>OQ 

/ IV>j,+ (z) 

J r 



■ dx 



< +oo, j,k = l,2. (5.81) 



Indeed, we will prove next that (5.81) follows from (5.80). For simplicity we consider 
the case j = k = 1. (The proof for the remaining combinations of indices j, k 
proceeds analogously, cf. also [4]). Assuming (5.80) for j = k = 1, there exists a 
constant C > such that 



f- f-X f-OO 

/ <Pi,+ ( x ) / ipi,-(y)f(y)dy + ipi-(x) / ipi,+(y)f(y)dy 

JR J -oo Jx 

<C [ \f(x)\ 2 dx. 

JR 



dx 



(5.82) 



For fixed r G R and / G L 2 (R) satisfying f(x) = 0, for all x > r, the inequality 
(5.82) implies (restricting the interval of integration) 



or 



|-0i i+ (x)| 2 dx 



f r ^,-(y)f(y)dy <cf \f(x)\ 2 dx. (5.83) 

J— oo J —oo 



Thus, choosing f(x) = tpi t -(x), for x < r and f(x) — otherwise (then clearly 
/ G L 2 (R)) one obtains (5.81) with j = k = 1. Since ||*±(x)|| 2 2 = \tpi,±(x)\ 2 + 
|V 2 ,±W| 2 , (5.81) yields (5.41). 

Finally, ||/C*(or)||c2 = ||c 2 proves the statement about the z- wave functions 

associated with q. □ 

The solutions ^±(z,x) in (5.38) are analogs of the familiar Weyl-Titchmarsh 
solutions in the context of self-adjoint Dirac-type operators (cf. [31, Ch. 3]). 
The following is a consequence of (5.38) and Remark 4.4. 
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Corollary 5.5. Assume q G Lj^R) and suppose the NLS_ case p = —q. Then 

a c (D(q)) D R, (5.84) 
a c (D(q)) D R, (5.85) 
a p (D(q)) (~l R = 0. (5.86) 



Proof. Relation (5.86) holds by Remark 4.4 which excludes the existence of an 
L 2 (R) 2 solution F of M(q)F = XF near ±oo for all A G M. To prove (5.84) we can 
restrict ourselves to the case in which p(D(q)) ^ 0. Pick Ao G R. Then Remark 
4.4, (5.38), and (5.86) imply that A i p{D{q)) U a p (D(q)). Since cr r (L>(q)) = 
by (5.28), it follows that Ao G a c (D(q)) by (5.14). Relation (5.85) is then obvious 
from (5.16). □ 

As a consequence of (5.84), our frequent assumption z G p(D{q)) (especially in 
the next Section 6) automatically implies z G C\R. 

Interesting restrictions on the permissible location of eigenvalues of J'-self-adjoint 
Dirac-type operators D(q) under strong additional constraints on q were recently 
derived in [25] . 

Remark 5.6. Given normalized Weyl-Titchmarsh-type solutions ^ ±(z , x , x ) of 
M(q) , i'(z,x) = z^>(z,x) for z G p(D(q)) satisfying 

1>i,±{z,x ,x ) = l, ze P (D(q)) (5.87) 

for some x G R, one can formally introduce associated Weyl-Titchmarsh m- 
functions as follows. Denote by E(z,x,x ) a normalized 2x2 fundamental system 
of solutions of 

M(q)^{z, x) = z*(z, z), zeC (5.88) 
at some xq G R, that is, S(z, x,xq) satisfies (5.88) for a.e. x G R and 

E(z,x ,x ) = h, zeC. (5.89) 
One then partitions E(z,x,x ) as 

E(z,x,x ) = (e(z,x,x ) Z(z,x,x o) ) = J$£$) , (5.90) 

where 9j(z,x,x ) and (j)j(z,x,x ), j — 1,2, are entire with respect to z G C and 
normalized according to (5.89). Then the normalized Weyl-Titchmarsh solutions 
if?±(z,x,Xo) can be expressed in terms of the basis (Q(z,x,x ) &(z,x,x )) as 

&±(z,x,xo) = Q(z,x,x Q ) +m±(z,x )$(z,x,x ), z G p(D(q)) (5.91) 

for some coefficients m±(z,Xo)- Clearly, m±(z,x ) are analytic on p(D{q)) and 
they are the obvious analogs of the half-line Weyl-Titchmarsh coefficients, familiar 
from {second-order scalar and first-order 2x2) self-adjoint differential and differ- 
ence operators (cf., e.g., [2, § VII. 1], [31, Chs. 2, 3]J. It is tempting to conjecture 
that appropriate boundary values ofm±(z,x ) as z approaches a(D(q)) encode the 
spectral information on D(q), but this is left for future investigations. 
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6. Transformation Operators for ^-Self-Adjoint 
Dirac-Type Operators 

The principal goal in this section is to construct transformation operators in 
L 2 (K) 2 that intertwine the jT-self-adjoint operators D(q) and D{qz^) corresponding 
to the Lax differential expressions M(q) and M(qi\ ) ) in the focusing NLS _-case 

and to use these transformation operators to relate the spectra oiD(q) and D(qi{'), 
the principal goal of this paper. 

In the following we always assume Hypothesis 5.3 and freely use the notation 
established in Section 5 for D(q) as the maximally defined J'-self-adjoint Dirac op- 
erator in the focusing NLS_ case and the Weyl-Titchmarsh-type solutions ^±(z, x), 
z G p(D(q)), established in Theorem 5.4. 

We start with an elementary but important observation. 

Remark 6.1. Since the two z -wave functions ^-(z) and ^+{z) of D(q) are linearly 
independent, 

W(H>-(z,x),V + (z,x))=c(z)^0, (6.1) 
all other (nontrivial) z -wave functions ^ satisfy 

= a*_(z) +/3*+(z) for some a, (3 G C\{0}. (6.2) 
In addition, as we will prove next, 

\\*(z,.)\\-}eL 2 (R). (6.3) 

Indeed, *(z) = a*_(z) + /3*+(z) and hence 

\Mz)\\£ < ||a|||* + (*)||<? - \(3m-(z)\\c*\~\ (6-4) 

(5.40), and the fact that ^±(z,-) G ACi oc (M.) 2 , yield the existence of constants 
C± > such that 

H^aOII^ = (||* ± (^x)|| c2 ||*(z,x)||- 2 1 )||vl/ ± (z,x)||- 2 1 

< C±||^±(^, ^H-, 1 , xgM. (6.5) 

By (5.39) this implies that all z-wave functions <J/(z) associated with q, except 
^±{ z )> satisfy (4.11). Hence, Hypothesis 5.3 guarantees the existence of z-wave 
functions ^(z) satisfying Hypothesis 4-2 for all z G p(D(q)). In particular, all but 
two z-wave functions of D{q) (viz., ^±(z)) satisfy Hypothesis 4-2 at z G p(D(q)). 

Without loss of generality we will restrict our considerations to the special case 
a = [3 = 1 in (6.2) for the remainder of this section up to (6.126), that is, we choose 

*(z) = V_(z) + V + (z) (6.6) 

in the following. 

Next, we pick some fixed z\ G p(D(q)). We take ^(z\) as in (6.6), ^(z\) — 
ty_(zi) + $f+(zi), where ^±(zi) satisfy (5.38)-(5.41) with z replaced by z\, and let 
JC^(zi) be the corresponding zT-wave function associated with q. By Theorem 5.4, 
K,^±(zi, •) satisfy (5.38)-(5.41) with z replaced by ~z~T. Moreover, 



fc*(*i)=£*_(zi)+K:* + (*i). 



(6.7) 
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Define the NLS potential q^ as in (2.38) and consider the zi-wave function 
^zi(zi) = (<f>^\zi),<f)^\zi)) T associated with as defined in (4.7), 

*£>(zi, *) =lm(z 1 )\\^(z 1 ,x)\\ c iJC^(z 1 ,x), (6.8) 
and the zT-wave function {z{) associated with qi\^ as defined in (4.8), 

a;) = -lm{ Zl )\\^{z u x)\\^{z u x). (6.9) 
We recall that according to (4.10) the new NLS_ potential is then given by 

q£>{x)=q{x)+A<fc\z 1 ,x)i }l {z 1 ,x). (6.10) 

Of course, both D{q) and D(c$} (associated with the differential expressions M(q) 

and M(qil^), respectively) are J'-self-adjoint by Theorem 3.5 and Corollary 2.9. 

In order to motivate the introduction of transformation operators one can argue 
as follows. Since 

T(z, *(*!), - -(i/2)(« - «i)/ 2 + ^(z^iz^, (6.11) 

one computes, for every z-wave function T(z, •) associated with q, 
T«(z, x) = r(z, a;, AC*(si))T(*, a;) 

= -{i/2)(z- z 1 )T(z,x) + $W(z 1 ,x)*(z 1 ,x) ± r(z,x) 
= -(i/2)(z - Z!)T(z,x) - ^(z 1 ,x)W(^{z 1 ,x),T(z,x)) 
= -(i/2)(z - Zl )T(z, x) - *W( Zl ,x)W{*- (z u x) + 9+{z!,x), T(z, x)) 

= -{i/2){z-z 1 )^T{z,x) + 2^{z l ,x) 



dx'*_(zi,a;') T criT(z,a;') - / dx' x') T cr 1 T(z, a;') 



(6.12) 

- $W(zi,a;)[W(*_(zi, -a), T(z, -a)) + a), T(z, a))], a > 0. 

To arrive at (6.12), we used the integrated form of (2.6). Replacing T(z) by F <E 
L 2 (WL) 2 in (6.12), noting that 



]hn.m£\W(y±(z 1 ,±a),F(±a))\ =0, Fe L 2 (R) 2 , 

afoo 



3.13) 



and repeating the same argument with T(z) replaced by fCT(z), then leads to the 
introduction of the following transformation operators T Zl and T— in i 2 (M) 2 , 

L 2 (R) 2 3 F(x) i ► (T Zl F)(a;) =F(a;) + 2$W(z 1 ,a;) (6.14) 

pX pOO 

x / dx'^-(z u x') T ^i F ( x ')- dx'^ + {z u x') T (j 1 F{x') , 

. J — oo J cc 

i 2 (R) 2 9 F(z) ~ (T-F)^) = Ffc) + 2K$V>(z 1 ,x) (6.15) 
x / dx' {K^-{z ll x')) T a 1 F(x l ) - / dx'(/C*+(z 1 ,a;')) T cr 1 F(x') , 

_J—oo J X 

zi e p(D(q)). 

That T Zl and T^- are in fact bounded operators in L 2 (K) 2 follows from Lemmas 
6.2 and 6.3 below. 
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Lemma 6.2 (Talenti [42], Tomaselli [44] (see also [4], [35]). 

Let f G L 2 {R), U G L 2 ((-oo, R]), V G L 2 ([i?,oo)) for all R G 

following assertions (i)-(iii) are equivalent: 

(i) There exists a finite constant C > such that 

dx U(x) dx'V(x')f(x') <C dx\f(x)\ 2 . 
JR. Jx Jr 

(ii) There exists a finite constant D > such that 



Then the 



(6.16) 



/ dx V(x) [ 
Jr J - 



V(x) / dx' U(x')f(x') 



<D [ dx\f(x)\ 2 . (6.17) 
Jr 



(Hi) 



sup 



dar|l/(a:)| : 



dx\V(x) 



< oo. 



Lemma 6.3. Assume Hypothesis 5.3 and Z\ G p(D(qj). Then the operators 

1 



dx' \\ i &±(z 1 ,x , )\\c>f(x'), 



dx' - 

T oo ||*±(«i,&')ll<? 



/(*') 



(6.18) 

(6.19) 
(6.20) 



(6.21) 



||*±(zi,a;)||c2 J± c 
L 2 (R)3 f(x)^\\y ± (z u x)\\ 

are bounded in L 2 (R). 

Proof. Using (5.38), (5.39), (5.41), and (5.67) one obtains 

sup(/ cfe||#_(zi,a:)||£2 J ( / dx \\^_( Zl ,x)\\c2 J < oo 

rER \J-oo J \Jr J 

and 

sup^y" dx\\^+( Zl ,x)\\-^ ^ dxW^+izux)^ <oo. (6.22) 

By Lemma 6.2, (6.21) is a necessary and sufficient condition for the operators asso- 
ciated with •) in (6.19) and (6.20) to be bounded in L 2 (R). Similarly, (6.22) 
is a necessary and sufficient condition for the operators associated with ^ + (z\, ■) 
in (6.19) and (6.20) to be bounded in L 2 (R). □ 

Thus, (6.19), (4.7), and (4.8) imply the following result. 

Corollary 6.4. Assume Hypothesis 5.3 and z\ G p(D(qj). Then the operators T Zl 
andT— defined in (6.14) and (6.15) are bounded operators in L 2 (R) 2 . 

Proof. By (6.8) one obtains \\<f> ( ^ {z u x) \\ £2 = \fm(zi)\\\*(zi,x)\\£. Applying (6.5) 
then yields the existence of constants C± > such that 



\^{z u x)\\ c2 = |Im(z 1 )|||*(z 1 , ^H-, 1 

< |Im(zi)|C , ±||*±(zi, a:)!!^ 1 , xeR. 



(6.23) 



At this point an application of Lemma 6.3 proves the boundedncss of T Zl in L 2 (R) 2 . 
Since ||/CC?||c 2 = l|G||c 2 for all G G C 2 , the same arguments prove boundedness of 
f w in L 2 (R) 2 . □ 
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In order to motivate the introduction of the inverse transformation operators, 
one inverts the matrix T(z, ^(zi), K,^{z\)) to obtain 

T(z, *(*!), /C*^))- 1 = c(z, Zl ) [ - (i/2)(z - zl)I 2 + (*i) x ] , (6.24) 

where 

c(z, zi) - -4(z - zx)- 1 ^ - zi)- 1 . (6.25) 
Thus, for z-wave functions T z ^(z,x) associated with q z \\ 
T(z, x) = T(z, x, /Ctt^!))- 1 ^ (z, x) 

= c(z, zi) [ - (i/2)(* - zT)T«(z, x) + AC*(zi, z)/C$£>(zi, x)^{z, x)] 
= c(z, zi) [ - (i/2)(z - zT)T«(z, x) - K*(zux)W{K*W(z u x), T«(z, a;))] 
= c(z,z 1 ){-(i/2)(z-zi)?%>(z,x) 

- (zi , x) + /C* + (z!, x)]^(/C$« (zi, a;), T« (z, a;)) } 

= 2*(z - z^-^Ti^z, a:) + 2AC*+(z 1 , a: ) da/ (/C$«(z 1; a^^iT^s, x') 

- 2/C*_ (zi , or) J dx' (/C$ « (z! , a/)) T ( r 1 T« (z, a;') 

- c(z, zi) [/C*_(zi, a ; )W(/C$«(z 1 , -a), T«(z, -a)) 

+ /C* + (zi,a:)^(/C$W(zi, a ),T«(z,a))], a > 0. (6.26) 

Replacing Tfi^(z) by F e L 2 (R) 2 in (6.26), and repeating the same argument for 
JCT z \\z) instead of T z \\z) , then leads to the introduction of the following (inverse) 
transformation operators S-^ and S Zl in L 2 (R) 2 , 

L 2 (R) 2 3 F(x) ^(S^F^x) = F(x) 

/X 
dx' (lC^(z u x')) axFix 1 ) 
-oo 

/■OO 

-2/C*_(zi,.t) / da;' (/C$W(zi,a/)) (TiF(a:'), (6-27) 

J a; 

L 2 (K) 2 9F(x)^(S Zl F)(a;)=F(x) 

+ 2*+(z 1 ,a;) j dx'^{z 1 ,x') T a 1 F{x') (6.28) 

pOO 

-2*-(z lt x) / dz' a^ifV), «i € p(D(g)). 

J a: 

Equations (6.20), (4.7), and (4.8) then imply the following result. 

Corollary 6.5. Assume Hypothesis 5.3 and z\ 6 p(D(qj). Then the operators S— 
and S Zl defined in (6.27) and (6.28) are bounded operators in L 2 (M.) 2 . 

Proof. Using again the estimate (6.23), one can follow the arguments in the proof 
of Corollary 6.4. □ 

We define two vectors F,G e L 2 (R) 2 to be J'-orthogonal if 

(F,JG) L 2=0, (6.29) 
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and we then write 

FLjG. (6.30) 

Since 

(F(x), JG(x)) C 2 = (JG{x))*F{x) = G(x) T (nF(x), (6.31) 

F± jG is equivalent to (F, JG) L 2 = [ dx G(x) 1 ' a 1 F{x) = 0, F, G G L 2 (K) 2 . 

Jr 

We also introduce the following notation of the ^7-orthogonal complement V ±J to 
a subset VcL 2 (M) 2 , 

V ±J ={Fe L 2 (R) 2 | FLjG for all G e V}. (6.32) 

In analogy to the orthogonality property of eigenvectors corresponding to different 
(necessarily real) eigenvalues of a symmetric operator in some complex Hilbert 
space H, one infers the ^-orthogonality of eigenvectors corresponding to different 
eigenvalues of a ^-symmetric operator S in H. Indeed, Sfj — zjfj, j = 1, 2, with 
z\ ^ Z2 implies 

MJh,h)n = {Jh,Sh) n = (JSfi,f 2 )n = (Jzifi, h)n = MJfu h)n 

(6.33) 

and hence 

(zi - Z2)ih,Jh)u = 0, implying {f u Jf 2 )n = 0. (6.34) 

Next, let do be an isolated subset of a(D[q^y) in the sense that (To can be 
surrounded by a positively oriented, rectifiable, simple, closed countour j ao C 
p(D(qi\^)) separating cr from the remaining spectrum <t(.D(<7z^))\<to. Then the 
Riesz projection onto the spectral subspace corresponding to cr is given by 

P ^ {(To) = "2^ / dZ Zh ^- (6 ' 35) 

The spectral subspace E^(ao) corresponding to cr is then the range of the Riesz 

projection, E^Vo) = M^W- 
We record the following result. 

Lemma 6.6. Assume Hypothesis 5.3. 

(i) Let zq,Sq G cfp{^D{qzi)) , zq ^ zq. Then the L 2 (K) 2 -eigenfunctions <&i^(zo); 
^i(zo) corresponding to z and z , respectively, are J -orthogonal, 

^{zo)±j^{zo). (6.36) 

(ii) Let zq G <7 P (-D(g^)) and (To an isolated subset of a(D(qi\^)) which does 
not contain zq. Then the L 2 {R) 2 -eigenf unction ^zi(zo) corresponding to zq is 
J -orthogonal to the spectral subspace E^ (<7o) corresponding to a a , 

*W(zo)±jXW(*o). (6-37) 

Proof. Assertion (i) is clear from (6.33) and (6.34). To prove (ii), one chooses a 
positively oriented, rectifiable, simple, closed contour "f aa which separates ctq and 
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{zq}. For F € L 2 (M.) 2 one then computes 



L 2 



(pU(a )F,J*U(z )) L2 = -± [ dz (£>(«&)) zI 2 y 1 F,J^(z )) L 

= -^ij^ dz(F,j{D{ q V)-zI 2 )-^{z )) 

Here we used the fact that since D^q^) is .J-self-adjoint, so is (D(q^) — zh) ■ 

□ 

Remark 6.7. From Lemma 6.6 and Theorem 4-5 one concludes that 

^(z 1 )±jlC^(z 1 ) (6.38) 
since z\ is a nonreal eigenvalue. Thus, 

*%>{z 1 )e{lC*U(z 1 )} J -*. (6.39) 

Remark 6.8. By Corollary 6.5, the operator Szy, z\ G p(D(q)), is well-defined 
and bounded in L 2 (M.) 2 . However, for future considerations it is more appropriate 

to restrict it to the closed subspace {/C(f>i^(zi)} ±J7 . Hence from now on, we will 

denote by S— the restriction of S— to {/C$i^ (zi)} ±J , 

S * = S*\{ K *v (xi) } ± *> *i €/>(£>(?))■ (6-40) 
An elementary computation (based on (6.31)) then reveals that 
{S^G){x) = {S^G){x)=G{x)+2K^>{z 1 ,x) f dx' (IC^ ( Zl ,x')) T a^x'), 

J — CO 

Ge{OW( Zl )} ij . (6.41) 
Next, we prove several results leading up to the principal theorems of this section. 
Lemma 6.9. Assume Hypothesis 5.3 and z\ € p(D(q)). Then 

ran(T Zl ) C {lC$£>( Zl )} ±J , *W( Zl ) € ker(^), (6.42) 
(*£>(*!), J*£>(zi)) L2 -Im(z 1 )(2^(*+(z 1 ),*_(z 1 )))- 1 ^0, (6.43) 

T Z1 S^G = G- t £; J ?" { $}» *£>(*i), Ge{lC*<»{z 1 )} ± *, (6.44) 

S W T ZI F — F, F e L 2 (R) 2 . (6.45) 
Proof. For brevity, we introduce 

/x />oo 
dx'^>-(zi,x') T aiF(x')- dx' a/) T aiF(x'), (6.46) 

-co J X 

F e i 2 (M) 2 , x e K, 

such that 

(T Zl F)(a:)=F(x)+2*( 1 1 )(«i,a : )^(a:;F) ) FeL 2 (R) 2 . (6.47) 
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Using 

£ x (x; F) = ^{z 1 ,x) T cr 1 F{x) for a.e. x £ M, (6.48) 

one computes 

{T Zl F,JJC^( Zl )) L2 = (F,JK^(z 1 )) L2 +2{^(z 1 )^-,F),JIC^(z 1 )) L2 
= {F,JIC^(z 1 )) L2 -(Im(z 1 ))- 1 I dx{\\^{z u x)\\l 2 ) x Z{x-F) 

JR 

= (F,JJC<S>{\\z 1 )) L2 + (Imz 1 )- 1 f dx\\*?)( ZuX )f *( ZuX )T aiF ( x ) 

Jr. 

= (F.J^Wjy +MzT) [ dx\\^(z 1 ,x)\\lH(z 1 ,x)~ T o.Fix) 

JR 

= 0, FeL 2 (Rf. (6.49) 

To justify the integration by parts step in (6.49) one can argue as follows. By 
(6.23), the estimate 

< IMzx)! 2 !!*^^)^ 



dx'w^^x'^ywFix'^y 

+ / dx'W^+iz^x'^yWFix')^ 

J X 

< iMzx)! 2 !!*^^)!!-, 1 ^!!*-^^)!!-, 1 f dx' \yst_{z u J)\\c>W 
+ c + \\y + ( Zl ,x)\\- 2 1 / dx>\\y + ( Zl ,x')\y\\F(x')\y 

J x 



(6.50) 



yields 

\\^X^-)\\lM-F)\^L\^), (6.51) 
using (6.3) and (6.19). Thus, 

Yim^d\\^{z u x)\\ 2 c2 Ux-,F)\=Q (6.52) 

(actually, lim^ioo • • • | — in (6.52) since all Lebesgue integrals involved in (6.49) 
are finite), which was to be proven. Hence, one concludes ICQ 1 ^ (zi) _lj-ran(T Zl ) . 
Next, using (2.6), Lemma 2.2 (iii), (6.8), (6.39), and (6.41), one computes 

(S w ^( Zl ))(x) = *W{z u x) + 2£*(*i,a0 f dx'(/C$W(z 1 , 2 ;')) T ( T 1 $( 1 )(z 1 , a ;') 

J — oo 

= $ « (Z! , z) - 2/C* (z! , a:) (2Im(zx ) ) - 1 W (/C$ « , z) , $ ( ^1,*)) 
= $ W (Z! , x) - /C* (zx , a:) (Im(z! ) ) - 1 1 1 $ W (zx , x) 1 1 ' 2 
- ^(z^x) - Im(z 1 )||vl/(z 1 , a ;)||- 2 2 /CvI/(z 1 ,a ; ) 

= (6.53) 

and hence G ker(5^ r ). 

For the proof of (6.43) we next assume G 6 {/C'&i^Z!)}^" 7 . Using 

2Im(z)* T (zi,a;) T cri/C*(zi,a;) = W(^ T {z!, x), JC^(z, x)) x (6.54) 
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(cf. (2.6)), one then computes 

(T Zl S-G)(x) = (S-G)(x) + 2$«(z 1 , 2 ;)£(x;S-G) 

/X 
dx' (lC$i) ) (z 1 ,x')) T a 1 G(x') + 2$ < z 1 1 ) (z 1 ,x)t(x;G) 
-OO 

+ 2(Im(z 1 ))- 1 $W(z 1 , a; ) (6.55) 
x f dx'W{^-{z u x'),lC^{z\,x')) xl f dx" (K,^ i ) {z ll x")) T a 1 G(x") 

. J —OO J — OO 

/■oo rx' 

- dx , W(* + (z 1 ,x'),lC*(z 1 ,x')) x , dx" (K^{z u x")) aM*") 

•J x J — OO 

By (6.23) one estimates 

W{^(z u x),)C^(z 1 ,x)) f dx , (/C$i 1 1 ) (zl,a; , )) T ( T 1 G(a; , ) 

J Too 



< ||vM*i,z)IHI*Mllc- f dx'\\^(z u x')\\ C2 \\G(x')\\ c2 

J Too 



=poo 



< iim^oiGili^^!^)!!^!!!*.^^)!!^ + \\* + ( Zl ,x)\\ C 2] 

x f dx , \\* ± (z 1 ,x f )\\£\\G(x')\\ c >, 

J Too 



(6.56) 



and hence the left-hand side of (6.56) is in L 1 ((=Foo, R]) for all R G M by Lemma 
6.3. Thus, 



lim inf 

x^=poo 



W(^> T (z 1 ,x),)C^(z 1 ,x)) f dx' (K^(z 1 ,x')) T ( j 1 G{x ) 

J ^foo 



0. (6.57) 



An integration by parts, using (6.57) and 

(lm(z 1 ))- 1 ^(z 1 ,x)W(^(z 1 ,x),IC^{z 1 ,x)) = -K>&(z x ,x), (6.58) 
then yields 

(T Zl S-G)(x) = G{x) + 2*<£(z 1 , X )Z(x;G) - 2(Im(z 1 ))- 1 <i>« ( Zl , x) 



dx' W{$!-{z\, x'), W(z 1 ,x'))()C$£>(z 1 ,x'))a 1 G(x') 

dx' W(* + {z 1 ,x'),K^{zi,x')){K^{z l ,x')) T a 1 G{x') 
= G(x) + 2$£>(zi,x) 

dx'9_(z 1 ,a;') T criG(a;') - / dx' Q + { Zll x') T a^x') 

-oo J X 



(6.59) 



(Again, lim^ioo | • • • | = in (6.57) since all Lebesgue integrals involved are finite). 
Here we introduced the abbreviation 

e±( Zl ,x) = V±(z 1 ,x) + *{z 1 ,x)\\V(z 1 ,x)\\-?W{y±(z 1 ,x),lCy(z 1 ,x)). (6.60) 

Actually, using the Jacobi identity 

AW(B, C) + BW(C, A) + CW(A, B) = 0, A,B,CeC 2 , (6.61) 
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one infers 

*±(zi)T^(/C*(z 1 ),*(z 1 )) + *(z 1 )^(*±(z 1 ),/C*(z 1 )) =£*(*i)W(*±(*i),*(*i)), 

(6.62) 

which in turn implies 

Q±(z u x) = (lm(z 1 ))- 1 W(*±(z 1 ,x),* T (z 1 ,x))<l>£>(z 1 ,x). (6.63) 
Combining (6.59) and (6.63) results in 

(T Zl S-G)(x) = G(x) + 2(lm(z 1 ))- 1 W^-(z 1 , X ), * + ( Zl , x))$<£> (zi, x) 

x f dx'^iz^x'y^Gix'). (6.64) 
Jm 

Thus, 

(T Zl S w G){x) =G(x) 

+ 2(Im(z 1 ))- 1 W(^-(z 1 ,x),^+(z 1 ,x))(G,J^ z \\z 1 )) L2 ^(z 1 ,x). (6.65) 

Applying (6.65) to G = ^(zi) G {lC$i\\z 1 )} ± - :r , one obtains the relation 

0=l + 2(lm(z 1 ))- 1 W^-(z 1 ,x)^ + (z 1 ,x))(^ i \z 1 ),J^(z 1 )) L2 , (6.66) 

which proves (6.43). Insertion of (6.66) into (6.65) proves (6.44). 

Finally, let F G L 2 (R) 2 . ThenT^F G (zi)} ±J by (6.42) and hence (6.41) 

yields 

{S^T Zl F){x) = F{x) + 2^{ Zl ,x)i(x-F) 

/X 
dx' ( Zl , x 1 )) a x [F(x') + 2* W (zi , (a/; F)] 

-OO 

= F(z) + 2$«(z 1 ,a;)£(z; J F) +2KM(z 1 ,x) f dx' ( Zl ,x')) T a^ix') 

- 2(lm(z 1 ))- 1 JC^(z 1 ,x) f dx' (W(K*W {z lt x'), (z 1; x')))j(x', F), 

J —oo 

(6.67) 



using 



-2^) (z^)^^^ (6.68) 
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The estimate 

\\^(z 1 ,x)W(K:^(z 1 ,x),^{z 1 ,x))ax; F)\\ c2 

<2|Im(z 1 )| 3 |jvl/(z 1 , a ;)||c 2 1 |e(x;^)| 
= 2|Im(* 1 )| 3 ||*(*i,z)llc* 



<2\hn{z 1 )\ s 



f dx>\\y_( Zl ,x>)\\ c 4F(x>)\\ c2 

J — CO 

p CO 

+ / dx'w^+izux'^ywFix'^y 

J X 

C-\\*-(zi,x)\\$ / ds'll^-^i.^llcHI^OHc" 

J — CO 

/>CO 

+ C + \\* + (z 1 ,x)\\- 2 1 / dx'\\* + (z u x')\\ C 2\\F(x')\\ C 2 

J X 



(6.69) 



then proves 

\\w( Zl )w(ic^Xzi),^X z iM--> F )\\ci e l2 W ( 6 - 7 °) 

by (6.19). An integration by parts in the last term ol (6.67), using (6.8), (6.9), and 
liminf \\K*(z u x)W(K.*<g{zi,x), *W(z u x))Z(x; F)\\ c2 

< ^Im^OllMJ^II^^i.ajJII^II^W^.arJll^l^F)! =0 (6.71) 

by (6.70) (in fact, lim^-oo | • • • | =0 in (6.71) since all Lebesgue integrals involved 
are finite), then proves S^T Zl F — F and hence (6.45). □ 



Next, we further restrict S— and define the operator S Zly — by 

S Zl ,— = S— | {$ (i) (zi)x<s >W (zi)} ± j = S ~\{^{ z ^),K^{z{)}^J- ( 6 ' 72 ) 

Lemma 6.10. Assume Hypothesis 5.3 and Z\ G p(D(q)). Then 

ker(T Zl ) = {0}, ran(T Zl ) = {$i 1 i )^ 1 ),/C$«(^)}^, (6.73) 

ker(5'-) = span{$« (*!)}, ran(S-) = i 2 (K) 2 . (6.74) 
Moreover, S Zl .— is the inverse ofT Zl , that is, 

ker(S Zl -) = {0}, ran(S Zli5r ) = L 2 (M) 2 , (6.75) 

S ZuzT T Zl =I 2 on i 2 (R) 2 , (6.76) 

T Z1 S ZUZI = I 2 on {^( Zl ),IC^(^)} ±S - ( 6 - 77 ) 

Proof. Suppose T Zl F = for some F G L 2 (R) 2 . Then (6.45) yields = S w T Zl F = 
F and hence ker(T Zl ) = {0}. The assertion ran(S l 2 I ) = L 2 (R) 2 in (6.74) is also clear 
from (6.45). Next, assume S— G = for some G G dom(S-) = {K,^{z{)}^ J (cf. 
(6.40)). Then (6.44) implies 

= T Z1 S^G = G - J°ff^ L * *<i>(*i), G G {K^Hzi)}^, 

(6.78) 

and hence G = c^i^(zi) for some c G C. This proves (6.75). 
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Next, suppose G € (zi), K,<S> { Z } {z-l)}^ . Then (G,J<S> { J?) L2 = and (6.44) 
imply T Zl S—G — G and hence 

{^{ Zl )X^{zi)}^ Q ran(T Zl ) C fr)}^ > ( 6 - 79 ) 

taking ran(T Zl ) C {/C$i^(^i)} J in (6.42) into account. The computation 

($<£( Zl ),JT Zl S^G) L2 = (*W( Zl ),JG) La (G,J^(zi)) L2 = (6.80) 

then proves (zi)-Ljvan(T Zl ) and 

ran(T Zl ) = {$W(z 1 ),/C$W(, 1 )} ±J (6.81) 

since ran(S w ) = L 2 (M) 2 . 

Finally, (6.76) and (6.77) are clear from (6.44), (6.45), (6.73), and (6.74). □ 

Next we state an auxiliary result. 

Lemma 6.11. Assume Hypothesis 3.2 and let £ e ACi oc (R), $ = ((f>i,(f> 2 ) T £ 
AC loc (R) 2 , F = (/i,/ 2 ) T edom(D{q)), and 

A_ e {G G ACioc(R) 2 | G,M(q)G e £ 2 oc ([-oo,oo)) 2 }, (6.82) 

A + e {Ge AC loc (M) 2 |G,M( (Z )GeL 2 oc ((-oo,oo]) 2 }- (6-83) 

TTien, 

(M(g)f$)(x) =ie , (x)(0 1 (a;),-0 2 (a;)) T + f (a;)(M(g)*)(x) /or o.e. i£l (6.84) 
and 

/" dx'A_{x') T a 1 {M(q)F)(x')=r dx' (M(q)A_){x') T cnF^') 

J — oo J —oo 

+ iA_(x) ± F(x), (6.85) 

/■OO /*OC 

/ dx' A + (x') T a 1 (M(q)F)(x') ^ dx' (M(q)A+)(x') T a^ix') 

J X J X 

-iA+{x)^F{x) (6.86) 

/or a// a; G R. 

Proof. Assertion (6.84) follows directly from the definition of M(q). To prove as- 
sertion (6.85) one integrates by parts and obtains 

f dx'A_(x') T a 1 (M(q)F)(x') = i[a 2 ,.(x)f 1 (x) - a h .(x)f 2 (x)} 

J — OO 

+ i f dx' (-oi -qfi - a 2 -qf2- + a[_f 2 - a' 2 _fi){x') 

J —OO 

= f dx' ((M(q)A_)(x')) T a 1 F(x')+iA_(x)- L F(x), (6.87) 

J — oo 

using liminf^i^oc \A-(x) ± F(x)\ = (actually, lim^j-oo | • • • | = 0). Relation (6.86) 
is proved analogously. □ 

The following result shows that T Zl and Szt intertwine D(q z ^) and D(q). 
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Lemma 6.12. Assume Hypothesis 5.3 and z\ G p{D{q)). Then, 

D(qW)T Zl =T Xl D(q), (6.88) 

S^DfaM) = D(q)S w . (6.89) 

Proof. Using formulas (4.10) one infers 

M(q^)-M { q) = 4i( ° (6 . 90) 

Relation (6.84) applied to M(qW) with $ = ^(zi) and (cf. (6.46)) 

/X pOO 
dx'*_(z 1 ,a; , ) T criF(a;') - / dx' ^ + (z 1 ,x') T a 1 F(x'), 
-oo J X 

Fedom(D(q)) (6.91) 

then yields 

(M(<&>) (£(•; F) $W ( Zl ))) (*) - J *(z 1 ,x) T a 1 F(x)(^ i (z 1; x), -^ l (z 1 ,x)) T 

+ z 1 $W(z 1 , a; )e(x;F) (6.92) 

since M^)*^^) = Zi^^i). By (6.85), with A_ = #_(zi) and A+ = 
\I/_l_(zi), one infers 

px r-x 

I dx'V±{z 1 ,x') T a 1 (M(q)F){x') = zi / da; / tf±(;8i,a;') T CTiF(a;') 

«/ ±00 ±oo 

+ i*±(zi,a;) ± F(a;), F E dom(D{q)) (6.93) 
since (M(g)*±(zi)) T = z 1 ^ ± (z 1 ) T . Thus, 

(M(^))T Z1 - T Zl M(q))F = (M(g£>) - M(g))F 

+ 2M( <z W)($( 1 )(z 1 K(-;F))-2<i>«(ziK(-;M( (Z )F) 

- 4z( - ^« i (z 1 )^(z 1 )/ 2 ,0W i (z 1 )^(zi)/i) T 

+ 2* [tt(*i) T <riF] (z0, (z x )) T + 2z 1 $« (z^-; F) 

-2z 1 *<»(z 1 )£(--,F)-2i[*(z 1 ) ± F] {<fclM)^lM)) T 
= 0, Fe dom(D(gr)). (6.94) 

This computation also proves that 

T Zl dom{D(q)) C dom(.D(qrW)) and hence T Zl D(q) C £>(g£))T Zl . (6.95) 

Next, choosing G £ dom(L>(g)S*^-) (i.e., G G £ 2 (K) 2 such that S^G G dom(D(gr))), 
(6.95) implies 

S w D(q^)T Zl S^G = S w D(q£)[G- ($£>(*i), J^)~ L l(G, J*$(zi)) L3 ] 
= S-D{q^)G = S-T Zl D(q)S-G = D(q)S-G, Gedom(D(q)S-). (6.96) 

Here we successively used (6.44), (z x ) = (cf. (6.42)), and (6.45). Thus, 

one concludes 

dom(£>(g£>)) C dom(D(q)S zT ) and hence S w D(q^) C Dfe)^. (6.97) 
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Hence, 

dom(£>( g £))T ai ) = {Fe L 2 (R) 2 \ T Z1 F E dom(£>(<£))) } C dom(D(g)) (6.98) 

since S^(T Z1 F) = F E dom(D(qj) for F E dom(D(q2?)T zl ) by (6.97). Combining 
(6.95) and (6.98) then proves (6.88). Equations (6.88), (6.44), and (6.45) in turn 
imply 

S-D(q^)T Zl S-= S-D{q^)[l - (^{z^ 

= S^DiqW) = SjhT^DMSzt = D(q)S zT (6.99) 
and hence (6.89). □ 
Remark 6.13. One can prove, similarly to Lemma 6.12, that 
D{q^)f- = f-D(q) and kcr(T-) = {0}, 

S Xl D{q$) = D(q)S Zl and ker(S Zl ) = span {/C$«(^)} (6.100) 
and that TW- and S z I < ,,, 1 ±~ are transformation operators (inverse to 

each other). 

Summarizing the results of Corollaries 6.4 and 6.5 and Lemmas 6.9-6.12 obtained 
thus far, we are now in position to state one of the principal results of this section. 

Theorem 6.14. Assume Hypothesis 5.3 and Z\ E p(D(q)). Then T Zl and 
are bounded linear operators in L 2 (R) 2 which intertwine the operators D{q) and 

D{q%), 

D(qW)T Zl =T Zl D(q), (6.101) 

SzrD(qV>) = D(q)S w . (6.102) 

Moreover, 

kcr(T Zl ) = {0}, ian{T x J = {zi)} ±s , (6.103) 

kcr(S-) = span {*£>(*!)}, ran(S^) = L 2 (M) 2 , (6.104) 

and S Zl _— (cf. (6.72)) is the inverse ofT Zl , that is, 

S Zl .-T Zl =I 2 on i 2 (R) 2 , (6.105) 

T Z1 S Z1 - = I 2 on {^(zi),^^)}^. (6.106) 

An analogous statement can be formulated when z\ is replaced by zj, that is, for 
the operators T— and S Zl , but we chose not to dwell on it here due to the symmetry 
of the arguments. 

We conclude with the principal spectral theoretic result of this paper. 
Theorem 6.15. Assume Hypothesis 5.3 and 2 Z\ E p{D{q)). Then 

<j(D(qV)) = a(D(q)) U {z u z^}, (6.107) 
ffp (£>(«£>)) = a p (D(q)) U {z u zi}, (6.108) 
a e {D( q W)) = a e (D(q)). (6.109) 



2 We recall that z e p{D(q)) implies z G C\M by Corollary 5.5. 
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In other words, constructing the new NLS_ potential qi\ amounts to inserting a 
pair of complex conjugate (nonreal) L 2 (R) 2 -eigenvalues, z\ andzj, into the spec- 
trum of the background operator D(q), leaving the rest of its spectrum invariant. 

Proof. We will prove (6.108) and (6.109) from which (6.107) follows by (5.27). By 
Theorem 4.5 one has zi,zi <G cr p (Z)(gr|^)) . Thus, (6.108) is equivalent to 

a p (D( q ^))\{z 1 - 1 } = a p (D(q)). (6.110) 

Next we denote 

X( 1 ) = {$W(z 1 ),/C$W(z 1 )} ±J . (6.111) 

Since D{q^) is J'-self- adjoint, X^ is a closed, D(qi\^) -invariant subspace. In 
addition, we denote by D^q^) | (1) the part of D{q^) in with 

dom(£>(g«)| x(1) ) =X«ndom(D( g «)). (6.112) 

From (6.101) and Theorem 6.14 it follows that 

D(qW)\ xll) = TD(q)T^ 1 on X™ (6.113) 

and 

D(q) = T~ 1 D(q'^)\ x(1) T on L 2 (R) 2 , (6.114) 

with T = T 7 and T _1 = SWl , . . , . ,1, (cf. Lemma 6.10). Assuming 

I *ki (»i),/C*ki («i) / 

[i e a p {D{q)) and 

D{q)F{ l i) = nF(fi), + F(p) e dom(£>(g)), (6.115) 
then (6.113) implies 

i?(gW)| JC(1) TF(/i) = /iTF(/i). (6.116) 
Since kcr(T) = {0} (cf. (6.103)), this implies 

<Tp{D{q)) C <7 p (D( g «)) and hence tr p (I>(«)) C a p (£>( g «))\{^,ir}. (6.117) 
Conversely, assuming v e cr p (D(gi^))\{zi, zj} and 

B^^w^M^^'M. O^FWMedom^^))), (6.118) 
one infers from (6.114) that 

D{q)T~ 1 F^{u) = vT^F^Xv). (6.119) 
Since kcr(T- 1 ) = {0} (cf. (6.75)), one concludes T" 1 ^ 1 )^) ^ and hence 

a p (D(q))Da p (D( q ^))\{ Zll zT}, (6.120) 

implying (6.108). 

Next, one observes that (6.114) also implies 

a e (D(q))=a e {D(q^)\ x(1) ). (6.121) 

To prove (6.109) we note that (6.121) implies 

a c (D(q)) = <T e (D(qW)\ xW ) C a e (D(q^)). (6.122) 

Conversely, let A e cr c (D ) ) and {G„} n eN be a singular sequence of D(qi\^ 
corresponding to A, that is, a bounded sequence in dom^D(qi]^)) without any 
convergent subsequence such that lim„^oo (Z)(gi^) — Xl)G„ = 0. Since L 2 (M) 2 
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permits the direct sum decomposition (not to be confused with an orthogonal direct 
sum decomposition) 

L 2 {R) 2 = X^+span{^(z 1 ),)C^(z 1 )} (6.123) 
(here one uses the fact that $i^(zi)_Lj/C ( f>i^(zi), cf. (6.38)), one can write 

(^( Zl ),J^( Zl )) L2 (K^{ Zl ),JK,^{ Zl )) L2 

(6.124) 

where F n G n dom(Z)(g^)) is the projection of G n onto the space 

the projection being parallel to the subspace span {^i^ (zi), K,&i\\ Z i)} . Since the 

coefficients of $zi(zi) and fC<S>i\ (zi) in (6.124) are bounded with respect to n G N, 
one can assume that they are convergent by restricting to a subsequence. Thus, for 
some real numbers c\ and C2, one has 

(D(q^)-M)F n - Cl ^( Zl ) + c 2 IC^( Zl ). (6.125) 

Since A« is Z) (^-invariant, F n G X« implies (Dfa^) - A/)F„ G X«, and 
since is closed, the limit Ci^zi(zi) + c 2 /C$i^(zi) also belongs to X«. But 
this implies c\ = c 2 = 0, because (zi), JX$i^ (zi)) L2 = (cf. (6.38)) and 
(^(a*), JS^i)^ ^ 0, (/C^^i), J-AC^i^zi))^ ^ (cf. (6.43)). Thus, one 
concludes that {-F n }neN C X^ is a singular sequence of D(c$} corresponding 
to A since it is bounded and has no convergent subsequence. (Otherwise, by (6.124), 
{<j n } nG N would have a convergent subsequence, contradicting the assumption that 
it is a singular sequence). It follows that A G (T c (D(gij' ) ) \ X{1) ) = &e{D(q)) and 
hence a e (D(q)) = a e (D(qi^)\ x(1) ) D (T e (D(q£>)), proving (6.109). □ 

Theorems 6.14 and 6.15 are new. We note that they are proven under the optimal 
assumption q G L 1 1 oc (IR) (but they seem to be new under virtually any assumptions 
on q). 

In the special periodic case where the machinery of Floquet theory can be ap- 
plied, the issue of isospectral Darboux transformations is briefly mentioned in [33, 
Theorem 3]. This excludes the insertion of eigenvalues as in Theorem 6.15. Insert- 
ing eigenvalues into the spectrum of a self-adjoint one-dimensional Dirac operator 
(not applicable in the present NLS_ context) was investigated by means of trans- 
formation operators (along the lines of [18]) in [43]. 

We conclude this section with a few facts on iV-soliton NLS_ potentials. As 
shown in Lemma 4.6, the insertion of pairs of complex conjugate eigenvalues into 
the spectrum of the background operator D(q) can be iterated. To fix the proper 
notation, we now slightly extend our approach of Section 6 and we consider a 
more general linear combination 1 t> lk ( Z ,x) of &±( Z ,x): Assuming Zk G C\R, k = 
1, . . . , N, one defines 

% h (z k ,x) = *-(z k ,x)+'y k * + (z k ,x), 7fe GC\{0}, k = l...,N (6.126) 

(as opposed to to our choice 7 = 1 in (6.6)). In obvious notation one then de- 
notes the corresponding Nth iteration of the construction of qi^ix) presented 
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in Sections 2 and 4 (cf. (2.38), (4.10)), identifying V(zi,x) and ^ 7l (zi,x), by 

Qzi,}..ZN,Jl,---,JN i X )- 

In order to describe a well-known explicit formula for gif^..z JVl7ll ... l7JV (x) (cf., 
e.g., [34, Sect. 4.2], [36], [41]) one introduces the quantities 



(6.127) 



and the 2N x 2N Vandcrmonde-type matrices 



V 2N (x) 



(I Zl 

1 z N 
1 z N+1 



N-l 

z l 



N-l 



-N 



<Pl 



Zlipi 



zn^Pn 



Z2N 



and 



V 2N (x) 



Vi 

/l z 1 .. 

1 z N 

1 Z N+i . . 



N-l 

Z N+1 ^Pn+i Z N+ iip N+1 



Z 2N 1 V 2N Z2N^2N 



z n 1( Pn 



N-l 

z n+1 <Pn+i 



JV-1 



JV-1 



"AT 



Pi 



z N ip N 



\1 2 2 Af 

The iVth iteration 



iV-1 

z jV+l Z^+lfN+l 



,N-1 



°2N P2Ar Z2Nf2N 

z w ,7i,...,7N (a;) is then explicitly given by 



JV-l 



y N 
N 



7 N 
Z N+1 



Z N I 

Z 2N ' 



f2N ) 



(6.128) 



(6.129) 



IZ!,... 



,7JY 



(x) = q(x) — 2i 



.det(y 2N (x)) 



(6.130) 



Denoting by D(q, 



det(V 2N (x))' 

i^.' ) .. ZN ,7i,...,7 N ) the associated jT'-self-adjoint operator in L 2 (K) 2 , 
repeated application of Theorem 6.15 then yields 

^(<?£L Wl7l ,.., 7 J) = a(D(q)) U {z u zi, z N ,zW}, (6.131) 

*pP(</£L w>71 ,..., 7 J) = <r p (D{q)) U {zi,zl, . . . , zn, zn}, (6.132) 

°e(D(qW. ZN ^ u ... nN )) = a e (D(q)). (6.133) 
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